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PREFACE 


Introduction. Fractals are concrete objects — pictures — that students now see in 
books, on posters, and on computer, video, and television screens. Such pictures of fractals 
provide vivid illustrations of abstract mathematics, but they also require mathematics to 
be created and understood. Fortunately, the mathematics involved in creating and inves¬ 
tigating fractals coincides with, and can help in understanding, the theory that students 
majoring in the mathematical sciences need to master in the introductory courses in com¬ 
plex analysis, real analysis, and topology: iterations, sequences, convergence, boundedness, 
neighborhood, continuity, compactness, homeomorphisms, measures, and so forth. 

Purpose. To convince students of the close relationship between theory and practice, 
the present book demonstrates how the “abstract” topics taught in introductory courses 
in analysis and topology apply to “concrete” pictures of fractals. Besides presenting stu¬ 
dents with applications of mathematics at several levels, fractals offer several additional 
advantages: they are familiar to many different students; they demand no prerequisites 
from fields outside of mathematics; and they provide research projects with statements 
accessible to undergraduates. Reflecting the ubiquity of fractals, the present book may 
either serve as a main text in a seminar on fractals or on research for undergraduates, or it 
may serve as a supplement — as do the instructor’s handouts on special topics — in such 
courses as complex analysis, real analysis, and topology. 

Chapter 1 introduces the fractals called “Julia sets,” and develops plotting algo¬ 
rithms based upon complex variables. The material includes a new bound 
for the size of quadratic Julia sets and for Mandelbrot’s set, which involves 
only basic algebra and the triangle inequality with complex numbers. The 
new bound then leads to a justification of the Non-Attracting Fixed-Point 
Inverse Iteration Method that utilizes only elementary complex variables. 

Chapter 2 begins with a simple construction of von Koch’s fractal snowflake. 

The exposition then proceeds with an explanation of Hausdorff dimension, 
and with a new, elementary proof of the Hausdorff dimension of von Koch's 
snowflake, which utilizes only a counting argument and the basic topological 
concepts of compactness and of Cauchy sequences of continuous functions. 

To accommodate readers with various mathematical backgrounds, each chapter begins 
with examples at the level of a beginner in the subject. For intermediate readers, each 
chapter then reviews the relevant theory from the core of complex analysis, real analysis, 
or topology, as needed. Following the theory, exploratory term projects allow students 
to explore lesser known aspects of fractals, at a challenging level involving computing 
experiments and theoretical conjectures, in a way similar to Gleick’s account of Mitchell J. 
Feigenbaum’s discovery of Feigenbaum’s constant with an HP-65 pocket calculator. Fin all}’, 
each chapter proposes research problems particular to the topic under consideration, which 
give students opportunities to investigate practical and theoretical extensions of the text 
and of their courses; the proposed problems may lead to publishable original research. 
Also, each chapter contains extensive solutions to all the exercises — with explanations 
and intermediate steps — as well as hints and suggested strategies for the term projects. 


Computing. Throughout the text, exercises and term projects allow students to • 
experiment with “supercalculators, 11 which are pieces of software or hardware — for exam¬ 
ple, Hewlett-Packard Company’s HP-28C, HP-2SS, or HP-48SX hand-held calculators, or 
any computer equipped with adequate mathematical software — that link graphic, sym¬ 
bolic, and numerical capabilities to one another. Such experiments also demonstrate how 
mathematical theory and pictures of fractals reflect each other. 


Please note that the programs listed herein are only working 
prototypes, which appear here only to illustrate the ideas explained 
in the text and to document how certain displays of fractals were 
obtained, without any warranty of any kind. 


Acknowledgements. 1 thank my students in various courses at Eastern Washington 
University, in particular, Randy S. Kimmerly, Barbara Kuhl, and Holly A. Ross, and my 
colleagues, chemist Dr. John Douglas and topologist Dr. Benjamin T, Sims, for having 
studied from early manuscript versions of this book and for having tried the algorithms 
on several types of supercalculators, computers, and workstations. The preparation of the 
book then proceeded thanks to the dedication of Craig A. Hall at Eastern Washington 
University’s Computer Graphics Laboratory, who installed the necessary hardware and 
software while I was typing in the text. For the penultimate draft of the “TfcX script." 
I acknowledge the editorial assistance of my long-time friends Dr. Stephen P. Keeler, 
mathematician and manager at Boeing Computer Services, and Joyce de Vries Kehoe. 
writer. All errors remaining in the final draft are mine. The camera ready copy then took 
its final form thanks to the assistance of the staff of the Department of Mathematics at 
the University of Washington in Seattle. 


Cheney and Seattle, Washington, June 1990. 

Yves Nievergelt. 
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CHAPTER 1 



QUADRATIC JULIA SETS 
in Complex Analysis 


Summary. The present chapter reviews the theory necessary to define the concept of 
“Julia set" and to design an algorithm to plot Julia sets. The resulting pictures illustrate 
many concepts related to complex numbers, at levels from precalculus through complex 
analysis. 

Prerequisites. The treatment of filled Julia sets, in the second section, demands 
only a working knowledge of the arithmetic of complex numbers, as reviewed in the first 
section. The presentation of Julia sets, in the fourth section, requires a familiarity with 
the contents of a standard introduction to complex variables (such topics are reviewed in 
the third section). 
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Chapter 1 QUADRATIC JULIA SETS 


0. INTRODUCTION 


‘‘Julia sets" are examples of mathematical objects informally called ■‘fractals." Pic¬ 
tures of Julia sets, such as those in exhibit 1. are important, because they provide a 
graphic medium for communicating certain abstract but useful concepts from real and 
complex analysis. For instance. Julia sets illustrate the sensitivity to initial conditions of 
the convergence or divergence of some sequences of complex numbers. Such sensitivity 
forms a part of the informal concept of “chaos" and is sometimes called the "Butterfly 
Effect": weather is so sensitive to small perturbations that a butterfly flapping its wings in 
Tokyo may cause the weather later to change in Paris. This sensitivity renders long term 
weather forecasting practically impossible. For a historical account of the discovery of the 
Butterfly Effect by M.I.T.’s Edward Lorenz, see Gleick's Chaos (reference [38]). 
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Exhibit 1. 


(*) 

(a) - (d) Pictures of Julia sets, (e) Similarity under magnification. 


There are advanced mathematical surveys of the theory of Julia sets, for instance, 
those by Barnsley [28], Blanchard [29], Devanev [30], Devaney and Keen [32], and Falconer 
[33], [34], with abundant examples of images of Julia sets, notably in the works of Mandel¬ 
brot [35] and of Peitgen ei al. [36], [37]. In contrast, this chapter provides an introduction 
to Julia sets, beginning at the level of precalculus and ending at about the level of a rig¬ 
orous introduction to complex analysis. Because a treatise about Julia sets would require 
a substantial amount of theory, this chapter aims at only two goals: a definition and illus¬ 
trations of the concept of Julia set; and a rigorous — but simpler than otherwise available 
— development of an algorithm to plot Julia sets w r ith computers or supercalculators. 







§1 COMPLEX GEOMETRY AND ALGEGBRA 


3 


1. COMPLEX GEOMETRY AND ALGEGBRA 

The present section reviews those concepts from elementary complex analysis that 
form the mathematical basis of "filled Julia sets ," 1 at the level of precalculus and calculus 

1 . 1 . The algebra of complex numbers 

A definition in terms of the set R of all real numbers will demystify complex numbers. 

Definition 1. A complex number is an ordered pair of real numbers, often denoted 
by r = (j. y), with x 6 R representing the first coordinate of z. and y € R representing 
the second coordinate of z , as shown in figure 1 . The set of all complex numbers is called 
the complex plane and denoted by the symbol C (“blackboard" or "special roman" C j. 

Thus, the set of complex numbers is the ordinary two-dimensional real plane: C = R 

EXAMPLE 1 . (1,0), (0,1), (-3/,. \/2), and (jr,-e) are four complex numbers. 

Besides the representation of complex numbers with cartesian coordinates, there exists 
another representation, with polar coordinates. 



Figure 1 . Cartesian coordinates. (x,y), and polar coordinates, (r. 0). 

Definition 2 , The modulus (also called the magnitude) of a complex number 
z — (i,y) is the number |r|, also denoted by the letter r, and defined by 

r = |r| = -Jx 2 +V. 

The principal argument (or simply the argument) of a complex number 2 = (r.y) is 
the number Arg( 2 ), also denoted by B (the Greek letter “theta"), and defined by 

0 _ A ( ) _ / -Arccos( j/ v/j 2 + y 2 ) if y < 0, 

\ Arccos(x/-y/r J + y 2 ) if y > 0. 







4 


Chapter 1 QUADRATIC JULIA SETS 


Thus, the modulus r = |rj represents the distance from the complex number ; to the 
origin, and the argument 8 = Arg(-) represents the angle, in the range — ~ < $ < t. 
formed by the complex number z and the positive first coordinate axis, as in figure 1 . 

EXAMPLE 2 . 

( 2 . 1 ) The complex number z = (1,1) has modul us r = |( 1 , 1 )| = -y/l^ + l 2 - s/2 
and argument 8 = Arg(l. 1 } — Arccos(l/Vl 2 4- l 2 ) = Arccos( l/>/ 2 ) = r/ 4 . 

( 2 . 2 ) The complex number is(0,l) has modulus r = [{ 0 , 1 )| = \/ 0 2 + l 2 = 1 
and argument 8 ~ Arg(0,l) — Arccos(0/v / 0 2 + l 2 ) = Arccos(O) = r/2. 

(2.3) The complex number u> = ( 1 , — \/3) has modulus r = |(1, — %/3)f = v/1 ^ 3 = 

2 and argument 8 — Arg(l, — x/3) — ■ • ■ — -Arccos( V 2 ) — — ~/3, 

While the preceding definition and example proceeded from cartesian to polar coordi¬ 
nates. the converse operation converts polar coordinates back to cartesian coordinates by 
projecting a point z = (x.y) onto each of the coordinate axes: 

r = r'Cos(f?} and y s r • sin( 0 ). 

Thus, with rs |(i,y)| and 8 = Arg (x, y), as shown in figure 1 , 

(r,y) = r ■ (cos(^), sin(#)). 

EXAMPLE 3. For the complex number with modulus r = and principal argument 
8 = tt/ 4. we obtain x = r ■ cos(8) = y/2 * l/\/2 = 1 and y = r • sin(0) = y/2 ■ 1 /y/2 = 1 

The first advantage of complex numbers lies in their algebra, described in the following 
two definitions, which provides a concise notation for dealing with algebraic and geometric 
operations in the plane. 

Definition 3. The addition of complex numbers maps two complex numbers, u■ = 
(u.r) and 2 = (x.y), to their sum, denoted by w + z and defined by 

w + z = (u,v)+(x.y) = (u + x,r + y). 

Thus, the addition of complex numbers is the ordinary addition of points ("vectors") 
in the plane. Graphically, the sum ir + 2 = (u + x,r + y) lies at the fourth comer of the 
convex parallelogram with the other three corners at the points tr = (u.r), z = (x.y). 
and the origin, ( 0 , 0 ), as shown in figure 2 a. 

EXAMPLE 4. (1,2)+ (3,4) = (1 +2,3 + 4) = (3.7). 

REMARK 1 . (TRIANGLE INEQUALITY) Since the modulus of a complex num¬ 
ber equals the ordinary Euclidean distance from that number to the origin, the modulus 
satisfies the same triangle inequality as the distance does: for all complex 2 and u\ 


\z + U'| < |zj + M, 
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with equality ( |r + u'| = |r|4- |u'| ) if, but only if, there exists a non-negative real number 
L > 0 such that either z = Cu.' or u* = iz (with iz — f(j.y) - (fi.fy). as with 
the ordinary multiplication of a number by a vector in the plane). This means that 
\z + u'| = |r| + |u'| if. but only if. r and u’ lie in the same direction from the origin. 

Definition 4. The multiplication of complex numbers maps two complex numbers. 
u- = (u.r) and z = (x.y), to their product, denoted by tr ■ r. or tc x z (which does 
not denote a cross product in the present context), or irr. with 

w: - (u.v)(x,y) = (ux - vy, uy + t-j). 

EXAMPLE 5. (1, 2)(3,4) = ((Ix3)-{2x4),(lx4) + (2x3)) = (3-8,4 + 6) - (-5.10). 

Graphically, the multiplication of a complex number z by a complex number w 
rotates z by the argument of w and expands or contracts r by the modulus of u-. as 
in figure 2b. In mathematical symbols, 

Arg(irx) = Arg(tr) + Arg (c), 


Im 


w + z * (n+x,y+y) 




(b) 

Figure 2. Graphic interpretations of the complex addition and multiplication. 

A verification of the graphical interpretation of the complex multiplication consists of 
expressing the multiplication not with cartesian coordinates (as in the definition) but with 
polar coordinates (as in the graphical interpretation). To this end, let r = |r| and 
6 = Arg(r) denote the polar coordinates of z, and let s = |tr| and ( Arg(u') ( C is 
the Greek letter “zeta"') denote the polar coordinates of tr: 

x = r ■ cos(0) and y = r ■ sin(0), 

u = s*cos(() and r=5-sin((). 
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Consequently, 

w: - (u, v)(x. y) = (ur - vy, uy + rr) = 

(s ■ cos(0 ■ r ■ cos(0) - 5 • sinfO • r • sin(0) , s • cos(0 • r • sin(0) 4- s ■ sin(0 • r ■ cosi 9 j i 
= rs(cos(t,'} cos(0) — sin(C) sin(0) , cos( <,*) sin(0) + sin(() cos( 9 )) 

= rs(cos(C + 9 ), sinf; + B)). 

Thus, the product ivz = (u,r)(x,y) is the complex number with modulus |*i* - j = r$ = 
and argument Arg(u'z) =s ( + 8 =s Arg(u>) + Arg(r). 

EXAMPLE 6. With z — (y/3, 3) and u> s: (\/3,1) the product becomes 
zw = (\/3,3)(\Z3,1) = (v^v/3 - 3 xl,v / 3xl + 3xv^) = (0,4\/3). 


Verify that \zir\ = |r| • |u’{ and that Arg(zu’) = Arg(r) + Arg(u’): 
t-u’l ^ |{0.4v^)! = 4^ = 2^3 x 2 = |(^3,3)| • |(%/3,1)| = \z\ ■ |w|, / 

Arg(ru-) = Arg(0.4>/3) = tr/2 = 7r/3 + tt/6 = Arg(i/3, 3) + Arg (\/3.1) = Argfr) + 
Arg(u-). / 

With its addition and its multiplication, the set C forms a number field (also 
abbreviated as a field), which means that the operations of addition and multiplication 
satisfy the properties listed in table 1, the proof of which follows through straightforward 
algebra by expanding and comparing both sides of each equation. 


Table 1. The algebraic properties of the field of complex numbers. C. 

The following properties hold for all complex numbers (u, t 1 ), (x.y), and (p. q ). 


Associativity of + 
Commutativity of + 
Identity (0) for + 
Inverse for + 


((u.r) + (x,y)) + (p,?) = (u,r) + ((x.y) + (p.g)) 
(u.v)+(x,y)= (x,y)+(u.r) 

(x.y)+(0,0)=(x.y)= (0,0)+(x.y) 

(x.y) + (-x,-y) = (0,0) = (-x,-y) +(x.y) 


Associativity of x ((u,tt)(x, y))(p, g) = (u, t>)((x, y)(p, g)) 

Commutativity of x (u,u)(x,y) = (x,y)(u,r) 

Identity (1) for x (x,y)(l,0) = (x,y) = (l,0)(x,y) 

Inverse for x If (x,y) ^ (0,0) then (x,y)(x/(x 2 + y 2 ), -y/(x 2 + y 2 )) = (1,01 

Distributivity of x over + (u,t/)((x.y) + (p,g)) = ((u, v)(x, y)) + ((u, r)(p. q)) 


REMARK 2. (INCLUSION OF R INTO C ) The set of all complex numbers with 
second coordinate equal to zero, of the form (x,0), satisfies all the algebraic rules of the 
real numbers with ordinary addition and multiplication. Thus, 

(x,0) + (u,0) = (x + tt,0) and (x,0)(u,0) = (xti,0). 

Consequently, we often identify the set {(x,0) : x € R} C C with R. For instance, we 
identify (1,0) with 1 and (0,0) with 0. With such an identification, Ix(p.g) means 
(1.0)x(p. g) and, similarly, 0 + (p, q) means (0,0) + (p,g). 
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REMARK 3. (HISTORY) According to van der Waerden's A History of Algebra ,39\ 
pages 47-52, one of the first conceptions of complex numbers arose in attempts to solve 
polynomial equations in finance. At the time (in the fifteenth century}, complex numbers 
were considered to be “imaginary." before the availability of any formal definition. For 
such historical reasons, the second coordinate y of a complex number z — (x. y) is called 
the imaginary part (or the complex part) of abbreviated by y = Imfr). and the 
first coordinate, x, is called the real part of r, abbreviated by i = Re(r). Remember, 
however, that the "imaginary" part is also a "real" number, 

EXAMPLE 7. For z — (1.2), the real part is Ref 1,2) = 1 and the complex part is 
Im( 1,2) = 2. 

REMARK 4. (COMPLEX ARITHMETIC WITH SUPERCALCl'LATORS) The 
HP-28 and HP-4S supercalculators (as their predecessor, the HP-15C) can perform all 
the operations defined in the present subsection, as explained in table 2. 

Table 2. Complex arithmetic with the HP-28 and HP-48. 


Keys Comments Display 

(1.2) 1 ENTER I (3,4) 1 ± I Add two complex numbers. (4,6) 

(1.2) I ESTER I (3,4) l~x~ I Multiply two complex numbers. (-5,10) 

(1,2) f fih'TER I (3,4) r-Y1 Divide two complex numbers, (.44,. 08) 

(3,4) I ENTER 1 ABS Compute the modulus. 5 

(-1,0) I ENTER 1 ARG Compute the principal argument. 3.14159... 

(1,2) RE Extract the real part. 1 

(1,2) IM Extract the imaginary part, 2 

(3,4) C—*R Extract both coordinates. ^ 

5 [ ENTER ! 6 R—*C Build a complex number. (5,6) 


(l,l) R—*P (HP-28) Convert to polar form. (1.4142.7853...) 

(1,1) CEHIEE i POLAR (HP-48) Convert to polar form. (1.4142. 1. 7853...) 


REMARK 5. (LOCATION OF COMMANDS) On the HP-28C&S, the commands 
ABS, ARG, RE, IM, R—*C, C—*R, and R—‘P are in the CMPLX menu. 

On the HP-48SX, the commands ABS, ARG, RE, and IM are in the PARTS submenu 
of the MTH menu, the commands R-+C and C—»R are in the OBJ submenu of the PRG menu, 
and the command POLAR is available on the keyboard with the blue (right-hand) shift key. 

To get arguments in radians, set your supercalculator in Radians mode by means of 
the RAD command in the MODES menu. 
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Exercises 


Routine exercises 

Exercise 1. Perform the following operations (which illustrate the properties listed 
in table 1) either by hand or with a supercalculator. 


(1.1) (0.0) +(2.3) = 

(1.2) (2,3)+ (0.0) = 

(1.3) (1.0) x (2,3) = 

(1.4) (2,3) x (1,0) = 

(1.5) (2.3)+ (4,5) = 

(1.6) (4,5) + (2.3) = 

(1.7) (2,3) x (4,5) = 

(l.S) (4.5) x (2.3) = 

(1.9) (2,3) + (-2,-3) = 

(1.10) (-2,-3) +(2,3) = 

(1.11) (2.3) x (Vi3,- 3 /n) - 


(1.12) ( 2 /i3i - Vis) x (2.3) = ) 

(1.13) ((2,3)+(4,5))+(6.7) = 

(1.14) (2,3) + ({4, 5) + (6,7)) = 

(1.15) ((2. 3)(4, 5))(6, 7} = 

(1.16) (2.3)((4, 5)(6, 7)) = 

(1.17) (2,3)((4,5) + (6,7)) = 

(1.18) ((2, 3)(4, 5)) + ( (2, 3 )(6 ,7)) = 

(1.19) (-1,0) x (2.3) = 

(1.20) (2.3) x (-1.0) = 

(1.21) (1.0) x (1,0) = 

(1.22) (0,1) x (0,1) = 


The following exercises involve only straightforward computational proofs. 


Exercise 2. For each complex number z — (x,y) denote by 2 the complex number 
: = (x, —y), called the complex conjugate of z. Applying the definition of complex 
addition and multiplication, express the following three results in terms of x and y. 

(2.1) x + 2 = ( 9 ,?) (2.2) r — 2 = (?,?) (2.3) 2 x 2 = (?,?) 

Exercise 3. Prove that if a complex number z has modulus r and argument 8. 
then z 2 has modulus r 2 and argument 28. Thus, jx 2 ! = \z\ 2 and Arg( 2 2 ) = 2Arg(r), 


Exercise 4. Prove the reverse triangle inequality: for all u.' and : in C, 


ll*Mw||<k-w|- 


Exercise 5. For the purpose of this exercise, let 

. , , f 1 if y > 0, 
s.gnfy) = ( _j jf y < 0. 

Moreover, define a complex square root 

C —* C, z = (x,y) *—■ s/z = = (u,v) 

by the formulae [obtained by solving (x,y) = (u,t’) 2 = (u 2 - t ,2 ,2uv) for u and v ] 


u = 


r+ v / ?T 


and u = sign(y)< 


-x + yfx* + y 2 
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Prove that u’ 2 = z , that is, (u.v)(u,r) = (x.y). 

Exercise 6. For each complex number 2 £ C denote by it any square root of 
In other words, let tv represent any complex number such that w 2 = z. 

(6.1) Prove that — u* is also a square root of z: thus, ( —u - ) 2 = z. 

(6.2) Prove that there exists no square root of 2 other than u> and — u\ 

Exercise 7. Prove that if z £ C and if w =c V*. with the complex square root defined 
in exercise 5. then [tr| = y/\T\ (the ordinary real square root) and Arg(ir) = Arg(r)/2. 

Exercise 8. Either by hand or with the square-root key on a supercalculator, compute 
the square roots of the following complex numbers: 

(S.l) 2 = —1; (8.2) z = —4; (8.3) 2 = 1 = (0.1). 

Exercise 9. Prove the complex quadratic formula: for all complex coefficients (com¬ 
plex numbers) a, 6, c £ C with a ^ 0, all the solutions of the equation 

az 2 + bz + c — 0 


are 


— b + y/b 2 — 4a c 

- — - and 

2a 


— b — Vo 2 —4ac 
2a 


Exercise 10. Determine all the complex solutions of the equation z 2 + 62 4- 25 = 0. 


1.2. Sequences of complex numbers 

Many applications of complex numbers involve not only a few complex numbers but 
infinite sequences of complex numbers, usually denoted by (r* )*L 0 , or (u- ; , or simply 

(2*) for brevity. (Such notation as {2*} is also common for complex sequences, but we 
shall also need some notation for the singleton [set] consisting of only one number, 2*. and 
the only accepted notation for such a set is (rt); thus, the notation (r*) for sequences 
avoids a confusion with sets.) A sequence of complex numbers is a function $ : N —♦ C. 
with values usually denoted by z k instead of the functional notation s(k). The study 
of sequences of complex numbers begins with the question of whether a sequence remains 
bounded, diverges to infinity, converges to a limit, or none of these. 

Definition 5. A subset Z C C is bounded if, but only if, there exists a non-negative 
real number B such that 

1*1 < B 

for every element 2 £ Z . Of course, a subset Z C C is unbounded if, but only if. it is 
not bounded. Similarly, a sequence (2*) in C is bounded if. but only if, the set of all 
its elements, Z — { 2 * : fc € N}, is bounded. 

Definition 6. A sequence (2*) in C diverges to infinity if, but only if. 

lim 121 [ = 00. 

k — oc 
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Since the limit involves only real numbers, because |r* I € R->- . the limit means that for 
each real number B £ R there exists an index A's € N such that if k > Kg then 
|cjt| > B. Informally, the definition means that we may arrange that all moduli jz* i be 
as large as we want, provided that we exclude the first Kg terms. 

REMARK 6. An unbounded sequence need not diverge to infinity. 

EXAMPLE S. Let z* = fc(l +{ — !)*) for each k £ N. Then 


Z*«fc(l + (-!)*) 


f 2k if k is even, 
l 0 if k is odd . 


Thus, (z* ) neither remains bounded nor diverges to infinity. 

Definition 7. A sequence (z*) in C converges for tends) to a limit z € C if. 

but only if, for each positive tolerance (number) t there exists an index K t such that if 
k > I\ t then jz* — z| < t. 

Informally, the definition means that we may arrange that all the terms z* be as 
close as we want to z, provided that we exclude the first A* t terms. 

Proposition 1. Every converging sequence remains bounded. 

Proof. Suppose that a sequence (z*) in C converges to a limit : £ C. The definition 
of convergence applied to t = 1 (or any other positive number) ensures the existence of 
an index Kf such that if A'i < k then |z* - z | < 1. For such indices k > A'i the 
triangle inequality shows that 


1**1 = \zk - Z + z\ < |Zk - -[ + N < 1 + l-t- 
For the other indices ( k < A’i ), of which there are only finitely many, let 

B 0 = max {[**| : k 6 {0. A'i}} . 

(The notation max means the maximum element in the set following max: thus. Bq 
represents the largest of |zo|, |*i|, ..Iz/Cjl-) Hence, |z*| < B 0 for every k < I\\. 
Finally, let 

B = max{Bo, 1 + |*|}, 

which ensures that |z*| < B for every term z* in the sequence (2*). □ 

REMARK 7. A bounded sequence need not converge to any limit. 

EXAMPLE 9, The sequence (z*) defined by z* = (-1)* remains bounded (it 
alternates between one and negative one) but it does not converge. 


2. FILLED JULIA SETS 


The ‘’filled Julia set" associated with a complex function determines whether an initial 
point 2 q € C gives rise to a bounded or an unbounded sequence (z*) through iteration 
of the function under consideration, as explained in detail below. To avoid unnecessary 
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technical difficulties and to focus on the nature of filled Julia sets, the present chapter 
restricts itself to those complex functions defined by quadratic polynomials, 

2.1. The concept of filled Julia set 

Notation. For each complex number c £ C let f c denote the function 

f c ■ C —* C. with /t(:) = ; 2 + c 


EXAMPLE 10, If c = 0 then / 0 (r) = z 2 + 0 = z 2 . Similarly, if c = t = (0.1) then 

/i(*) = z 2 t- 

Deflnition 8. For each natural number n € N. the n-th iteration of / c is the 
function defined inductively by 

/;■(.-) = / c ( 2 ) = .- 2 + c, 

!?(--) = /.(At--)) = h (.- ! + c) = (.- ! + c) 2 + c, 

/ t ° 3 ( = ) = U (/. (/«(«)» = (U ! + O’ + e)’ + e. 

/?"(») (/?"■■(*))• 

Next, for each complex number ro, define a sequence (:„) by = /°"(zo)- Thus. 

-) = /c( 2 o) = 2 o + c ' 

--3 = o) = /, U'M) = -•; + c = (r 0 2 + c) 2 + c, 

•"] = ‘"2 + C . 

*n-f 1 "f- C. 


By induction, the formula z n +i = z 2 + c specifies the entire sequence ( 2 „). Remember 
that each complex number z„ in the sequence depends not only upon n but also upon 
c and 2 q , but the notation omits such references for simplification. 

EXAMPLE 11. For c = -1 the polynomial / c takes the values f-i(z) = z 2 - 1. 
If 2 0 = 1 then the sequence (z„) becomes 1, 0, —1, 0, —1, 0, —1, ..., because 

20 = 1, 2 j = zl+c= 1 2 -1 =0, z 2 = :\ + c= 0 3 -l = -1, z 3 = z| + c = (-1) 2 -1 =0. 

and so forth. Thus, with 2 0 = 1, the sequence ( 2 „) alternates between 0 and -1 
forever: in particular, the sequence (z n ) remajns bounded. By contrast, if zo = 2 then 

Zq = 2, 2 ] = Zq + c = 2 2 - I = 3. Z 2 = Z 2 + c = 3 2 - 1 = s, r 3 = + c = 8 3 - 1 = 63. 
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and so forth: with zo — 2, the sequence (~„) diverges to infinity (as proved below). 

Definition 9. The filled Julia set of f c is the set K c of all complex numbers 
r 0 for which the sequence (r„) = {f° n (:o)) remains bounded. 


EXAMPLE 12. For the particular value c = 0, the quadratic polynomial f c = / 0 is 
the squaring function, /o : C —♦ C, w r ith / 0 (r) = r 3 . Consequently, 


_ .2 

■‘' 0 ' 


2 1 =■ foi-o) — 
z 2 ~ /o(-1 ) — — ( 2 o) “ 

Z Z = /o ( 2 2 ) ~ 2 \ ~ t 2 *} 2 — : 




Figure 3, The sequence (r„) = (/° n (r 0 )) diverges to infinity if 
|zo| > 1 and it remains bounded if Jr 0 | < 1- 


Thus. \zn | — (ro) ,: f = i 2 ol <2 "'^ from which tw-o cases emerge: 


{ If |r 0 | > 1 then |z n | = (r- 0 ^ diverges to infinity. 

If \zq | < 1 then |z n | = Izol 12 "* < l (2n) = 1 remains bounded. 

The preceding considerations prove that the filled Julia set A'q of the squaring function 
/ 0 : C —> C, /o(z) = 2 2 consists of all points at a distance of at most one from the origin, 
as in figure 3: 


A'o = {r 0 : ro £ C and |z 0 ! < 1} . 
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The example of the squaring function, with c = 0. required only basic inequalities with 
moduli to determine its filled Julia sets: unfortunately, all values of c other than c = 0 
and c — —2 (examined in example 21 in the fourth section) require substantial analysis 
for determining the nature of the filled Julia set. Fortunately, a few partial results that 
involve only some algebra and inequalities provide some information about the location 
and size of the filled Julia set. For instance, the most elementary such result guarantees 
that each filled Julia set K c is bounded. 

Theorem 1. For each c £ C let R z = max{2.jc|) (the larger of 2 and jc j. 
Then the filled Julia set K c of the quadratic polynomial f c ij bounded by R c . -n the 
sense that if zq £ K c then |z 0 | < R c . Thus, no point in K c may have a magnitude 
greater than R e . 

Proof. (The present proof lacks elegance, but it may inspire you for exercise 12. at 
the end of this section, and for the project on Mandelbrot’s set, at the end of the chapter, i 
Suppose first that |c| < 2; in this case. R c = max{2. |cj} = 2. To show that K c is 

bounded, consider a point Zq £ C such that jzol > Rc and prove that z 0 gives rise 

to an unbounded sequence, as follows. Set d — |zo| — 2 > 0, so that [zq| = 2 + d. and 
estimate |r! |: 

1*11 = ko + c l - ikoi 2 - kl| = |(2 + d ) 7 - id) 

= |4 + 2 d+d 2 - |c|J = |2 + 2d + d 2 + (2- |c[)| > 2 +2d. 

Thus, if [z 0 | = 2 + d. then |rj| > 2 + 2d. By induction, the same calculation shows that 
|r n | > 2 + 2 n d, which diverges to infinity. 

If jc| > 2 then R c = max{2, jc|} = |c| but a similar argument leads to the same 

conclusion. If |zo| > R c = |c|, set d — |z 0 t — (c|, so that |zo! = |c| + d. Hence. 

1 * 1 1 = i *0 + e| > Ikoi 2 - |c|| = |(|c| + d) 2 - \c\\ 

= j|c| 2 + 2d|c| + d 2 - Jc|| = ||cj(jc[ - 1) + 2d\c\ + d 2 \ > |c|(2 - 1) + Ad = jc| + Ad. 

Thus, if |z 0 | — |c|+tf then |zi| > jcj + Ad. By induction, the same calculation shows that 
|r n | > |cf + 4 n d, which diverges to infinity. Consequently, if jz 0 | > Rc then zo £ K c . and. 
by contraposition, if z 0 £ K c then (rol £ Rc- D 

Notice that the theorem asserts only that if zo € R\ then |z 0 | 5 Rc\ it does not 
assert the converse, because there may exist points z 0 such that |2 0 i S. Rc hut z 0 £ I\ c : 
not all complex numbers with a magnitude of at most R c need be in the filled Julia set. 
In other words, the filled Julia set may be a proper subset of the disc with radius R c 
and center at the origin. 

REMARK 8. (A BETTER BOUND) The preceding theorem yields the standard 
estimate found in texts, R c = max{2, |c]} (see the chapter by Bodil Branner. "The 
Mandelbrot Set,” in the book edited by Devaney and Keen [32], page 80). Such an estimate 
need not be optimal, in the sense that the filled Julia set K e may be smaller than the 
bound R c . For instance, exercise 12 at the end of the present subsection shows that 

1 ± y/\ + 4 jc] 


2 
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provides a better bound for A’ c , because r c < and. still, if zo € A' c then |ro 1 < r f . 

The availability of a bound for filled quadratic Julia sets suggests a crude algorithm 
for determining whether a point zo € C lies in K c , 


Direct Iteration Method for filled quadratic Julia sets. For each 

c 6 C. the following algorithm generates points in the filled Julia set I\ c of 
the quadratic polynomial f c : C —* C with / c (z) = z 2 + c. 

Step 1. Choose a maximum number of iterations N € N, for instance. 
N = 12. depending upon time and accuracy (see Strang's text [40]. page 510). 

Step 2. Choose a bound for the filled Julia set A' e . for instance 

1 + v/1 4- 4jc| 

r < =-J-' 

Step 3. Color white all points zq € C such that |z 0 ] > r c : such points 
do not belong to A',-. 

Step 4. Test each point Zo such that |zo | < r c to determine whether 
zq 6 I\ c , as follows. Start computing elements of the sequence 

2i = 2 q + C, Z2 = Z-y + C, *3 — + C, . . . , r„ = Z n _ j + C, ... 

and stop if either of the following conditions holds. 

(oc) If at some stage n e the sequence returns an element 

such that | z n { > r c , then the sequence diverges to infinity (see exercise 
12.2); consequently, color the initial z 0 white. By symmetry (see exercise 
11), also color —zq white. Stop the iterations and select another zq- 
(A'e?) If n = N {the maximum allowed number of iterations) and if |z„l £ r c 
for all the computed elements r 0 ,..., zs then we do not yet know whether 
the sequence diverges or remains bounded; nevertheless, color the initial 
z 0 black (it may belong to the filled Julia set A' c ). Also color — Zo 
black. 


Exhibit 2 show's implementations of the Direct Iteration Algorithm on the HP-2S(C 
or S) and HP-48, with an example. 

Despite its simplicity, the algorithm just presented suffers from the disadvantage that 
while it can assert that an initial point z 0 gives rise to a diverging sequence, it cannot 
assert that a point zq belongs to A’ c . Thus, it colors all points Zq for which it cannot 
make a decision. Another algorithm, presented in the next section, provides a remedy for 
this disadvantage by producing only points in A' e . 
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Filled Julia sets Comments on the programs 

on the HP-28C&S for filled Julia sets 


Filled Julia sets 
on the HP-48SX 


« CLLCD 
FOR 

r- NEG r 


R -► C 1 

START 

1 2 

SQ c 


NEXT 


IF ft B S r 4 

THEN 

x w 

R -► C 

P I XEL 

END hi 


STEP hi 


STEP LC 
» 1 

&NTER 

1 scan 

■ STO 

« c FIBS 

* 

1 4" 1 

+ 2 

* r • STO r- 

32 / ■ hi 

* STO 

M 137 * 

2 / 

•s 1 STO 

s r R ■+ C 

PMftX 

s NEG 0 

R-eC 

PM I N 
» 

ENTER 

• setup 

* STO 

« setup 

0 r- 

sc*n & 

hi NEG 

R + C PM ft X 

s NEG r- 

hi -► 

NEG R-*C 

PM I N 

r- NEG hi 

NEG 

*c *n 
» 

ENTER 

•fill 

* STO 


TUTOR IflL 
.25 * c * STO 

USER fill 




Subroutine setup computes a 
bound for the filled Julia set, 

1 + y/l + 4|c| 
r = r c = - 2 -' 

and adjusts the screen so that it 
covers the square [—r,r] x'-r,r], 
hence also the filled Julia set; 
h equals the width of each pixel. 

ENTER and STOre in 'setup' 


For the HP-28 CSjS only, subrou¬ 
tine scan examines each pixel in 
either the upper or the lower half 
of the square [-r,r] x [—r, r]; scan 
and fill perform the same task 
as fill does on the HP-48SX. 

ENTER and STOre in 'scan' 


Program fill tests each pixel z 0 
= (j.y) in the square 

[—r. rj x [—r,r], by iterating 

M:) - : 2 + c 

twelve times within the loop 
START 

SQ c + 

NEXT 

If l*i 2 1 < r e then fill colors 
the pixels at r 0 and at —Zq in 
black. 

ENTER and STOre in 'fill' 
Tutorial: 

STOre a complex number in 'c' 
execute fill 
in the USER or VAR menu. 


* c ft B S 4 * 

1 ■+■ J* 1 2 

y- ' r- 1 STO r- 

64 / •h 1 STO 

ERASE r- 128 

y' 13 1 * r~ 

R -► C DUP NEG 

SWAP PDIM # 

131d # 128d 

PD I n 

» ENTER 

* * «• t up 1 STO 

« * e t ljf> 0 r- 

FOR w 
r- NEG r- 

F O R x x si 

R -► C 1 12 

START 
SQ c 

NEXT 

IF BBS r~ 4 

THEN x si 

R-*C DUP 
NEG PIXON 
P I XON 
END hi 

STEP M 

STEP 

* ENTER 
•fill * STO 

TUTORIftL 
. 25 * c 1 STO 

V ft R fill 



Exhibit 2 
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Exercises 

Exercise 11. Prove that every quadratic filled Julia set is symmetric with respect 
to the origin. 

Exercise 12. For each complex number c G C let 

1 + y/l + 4 I C I 

*- — _ -- — 

c 2 

(12.1) Prove that if |ro| > r c then (r n ) diverges to infinity. 

(12.2) Prove that regardless of |r 0 (, if there exists an element r* such that Jr*! > 
r f for some index Jt € N then (r„) diverges to infinity. 

(12.3) Prove that r c < max{2, |c|}. 

(12.4) Prove that if c is real and non-positive ( c €] — oc.O] ) then r c provides an 
optimal bound, in the sense that K c contains two diametrically opposite 
points each at distance r c from the origin. 

Exercise 13. Prove that ( 2 n ) = (/°"(*o)) cither remains bounded or diverges to 
infinity. (In contrast with other complex sequences, [for instance, the sequence in example 
8], the sequences (z„) that arise as iterations of the quadratic functions f c cannot be 
unbounded without diverging to infinity.) 


2.2. The Fixed-Point Inverse Iteration Method 


The Direct Iteration Algorithm demonstrated in the preceding subsection suffered 
from the disadvantage of not recognizing points in the filled Julia set. Another algorithm, 
demonstrated in the present subsection, generates points only in the filled Julia set. A 
further improvement of the new algorithm, explained in the next sections, will generate 
points only on the “edge" of the filled Julia set, which produces more reliable pictures. 
Both algorithms rely on the concept of “fixed pointT 

Definition 10. A fixed point of a function / : C —* C is a point c. such that 

/(-.) = ;.. 

EXAMPLE 13. Consider the function : C —♦ C with /j/*(c) = z 2 + V-i- To 
determine whether /]/ 4 has any fixed point, and. if so. to find such fixed points, solve 
the equation - r, which, in this particular example, yields 


/]/*( 2 ) — : - 

z 1 + = z, 

Z 2 - Z + 1/4 - 0. 


-(- 1 )± ■ 1 - 7 « 

2 ■ 1 


1 ± 0 


1 

V 


9 
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Hence, the function A /4 has one fixed point. 2 . = J /z. To verify this result, observe that 

/i/4U.)«/i/ 4 ( , A) = ( 1 /i) 2 + V< = 'A + J A = V: = thus. / 1/4 ( V2) = ‘A- ✓ 

EXAMPLE 14. For each complex number c ^ >/*, the quadratic polynomial 

/c : C —» C. / C (r) = z 2 + c 


has two distinct fixed points, 


1 + v / n r 4c , 1 - >/l - 4c 

=■>- 5 — “ d - 2 — 

To verify this assertion, substitute z + and z_ for z and check that / r (z + ) s= and 
fc( : ~) = z_, or apply the quadratic formula, which shows that 

if /«(*) = *, 

then z 2 + c - z, 

hence z 2 — z + c = 0, 

and 

_ —( —1) ± v'l 2 - 4 ■ 1 ■ c _ 1 ± ^1 - 4c 
— 2-1 ~~ 2 ’ 

Finally, since c ^ s / 4 . it follows that 1 - 4c ^ 0 and that the two fixed points 2 + and 
z_ are distinct ( z + ^ ), □ 

Every fixed point z, of a function f c lies in the filled Julia set A' c , because the 
sequence of iterations (z n ) starting at zq — z» remains at 2 , forever: Z] — /(z 0 1 s 
/(-.) = z.. *2 = /(zi) = /(z.) = 2 .. and so forth. r n = 2 , for every z n . Therefore. 
(z n ) remains bounded. Moreover, fixed points play a crucial role in the context of Julia 
sets, partly because they provide infinitely many points in the filled Julia set. 

Proposition 2. If 2 . u a fixed point of the quadratic polynomial f c then all 
its preimages under /*" (all the complex numbers z 0 for which there exists an integer 
n € N such that z„ = / e on (zo) — 2 . ) belong to the filled Julia set. 

Proof. Suppose that z„ = / c on (z 0 ) = z,. Since r, is a fixed point of f c . 

“» + l = fc(~n ) = /c(z» ) = Z,, Z n + 2 ~ fcl~ n+1 ) = /c(*• ) = Z., ... 

and all terms beyond z n equal z.. Thus, (z*) remains bounded, hence z 0 € A' e . □ 
The preceding proposition suggests an algorithm to plot filled Julia sets. 
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Fixed-Point Inverse Iteration Method. For each c£ C plot points 
in the filled Julia set as follows. 

Step J. Determine a fixed point of / c : choose either of the two fixed points 

1 ± VT~- 4c 
*■ = ** = ■ 

Step 2. Compute a preimage r _\ 

d: + c 

thus, chose either square root. 

z_! = ±y/z. ~ C. 


of the fixed point r. by solving 


and plot the point r_j. 

Step 3. Compute a preimage of the point just obtained by solving 

/e(*-a) = *-l i 


:1 2 + c = r_i, 


which gives 


r _ 2 = ±\A-i - c. 

and plot r_ 2 . 

Continue in this fashion until enough points appear on the plot. 


The experiments in the exercises reveal that all the preimages of a fixed point belong 
to the filled Julia set, but that they may fail to yield a satisfactory picture of the whole 
set: the preimages may stay confined in a subset so small that they do not oultine a repre¬ 
sentational picture of the filled Julia set. A slight variant of the inverse iteration method 
remedies this problem, but an explanation of that improved method requires concepts from 
topology, reviewed in the next section. For this reason, if you want to try a program at 
this point, then you may try the Non-Attracting Fixed-Point Inverse Iteration Method, as 
demonstrated in exhibit 4, page 33. in the fourth section. 
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EXAMPLE 15. Suppose that c = J / 4 . 

Step 1. Determine the fixed points by solving /j/ 4 ( 2 .) = r,. which gives 


from which it follows that /]/ 4 has only one fixed point, 2 . = l / 2 - Plot 
Step 2. Find a preimage of the fixed point by solving 

/ 1 / 4 U- 1 ) — z *-, 

‘ii + 1 /a - J /zi 

2-1 = ± y /~2 — ] A — = ± l h- 

Notice that one of the two solutions, J /j, coincides with the preceding point, 
2 ., Consequently, select the other solution, 2 _j = — V 2 * and plot z \- 

Step 3. Compute a preimage of 2 _) by solving 

fl/Ai *—2 ) — 2_i, 

* — 2 + V* = ~ V 2 . 

2-2 = ±\/- V 2 — _1 A - ±v /r TyvI = ±i\/3/2. 

Choose either solution randomly, for instance, 2_ 2 = tV3/2. and plot it. 

Step 4. Compute a preimage r _3 of for instance 


2-3 = -v/2_2-C = 3/2 — >/,, 

and plot it, and so forth. 


Exhibit 3 displays a picture generated 
by the Fixed-Point Inverse Iteration Method 
implemented on the HP-4SSX {see exhibit 4. 
page 33), Observe that all the points ap¬ 
peared to lie on the “edge" of the filled Julia 
set, a phenomenon explained in the fourth 
section. 
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Exercises 


Experiments and conjectures 

Exercise 14. The present exercise investigates the filled Julia set 

(14.1) Determine the two fixed points of /_; : C —* C, /_ 2 ( 2 ) = z 2 — 2 . 

(14.2) For each fixed point, compute the first three preimages, 2 _j, z_ 2 , and r_ 3 . 
as demonstrated in example 15. 

(14.3) Prove that if -2 < 2 _„ < 2 then -2 < r_ (n+1 j < 2 (with 
representing any solution of /_ 2 ( 2 _( n+ 1 j) s= 2 _„ ). Conclude that if the 
Fixed-Point Inverse Iteration Method starts from either fixed point of 
then all the points that it generates lie in the interval [— 2 , 2 ]. 

We shall prove later that the filled Julia set A ’_2 is indeed the closed interval [— 2 . 2 ], 
from —2 through 2 . However, the following two exercises illustrate how the choice of the 
fixed point may influence the quality of the picture produced by the Fixed-Point Inverse 
Iteration Method. These exercises also warn against quick unverified conjectures. 

Exercise 15. Let c = 0 and consider the squaring function, with /o(r) = z 2 ■ 

(15.1) Determine the two fixed points of the squaring function, /o- 

(15.2) For each fixed point of / 0 , compute the first few preimages as explained 
in example 15 with one modification: at each step, keep both preimages. 

— c and - \/£_n — c. Compare the plot resulting from one fixed point 
with the plot resulting from the other fixed point. 

Exercise 16. With c = — 3 /« let : C —* C, /_ 3 / 4 {*) — z 2 - ty»- 

(16.1) Verify that 2 . =s - J /2 is a fixed point of /_ 3/4 ■ and also determine the 
other fixed point. 

(16.2) Start from the fixed point, 2 , = — x h, and apply the Fixed-Point Inverse 
Iteration Method as explained in example 15, with one modification: select 
the positive square root at each step, so that 

‘0 = — Vji *-i = + %/— V2 - (— 3 A) — Vi' * - 2 = + \J i-i — (— 3 /»h 

To what limit 2 does the sequence of preimages (r_„) appear to converge? 

(16.3) Prove vour conjecture (prove that lim T ,_ 3 C ( 2 _„) = 2 ). 

(16.4) Start from the other fixed point and plot the resulting sequence 
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3. COMPLEX ANALYSIS AND TOPOLOGY 


3.1. The topology of the complex numbers 

The concept of limit applies not only to sequences, which axe functions from \ to 
C, but also to functions from C to C. 

Definition 11. Let / : D —+ C be a function defined on some subset D of the 
complex plane. Also, let zo € D represent a point in the domain, and let u'o € C. Then 
f(z) tends to Wq as z tends to z 0 if, but only if, for each positive tolerance f £ Ri. 
there exists a positive distance (number) d € R+ such that the following condition holds: 

If \: - z 0 \ < d, then z € D and |/(r) - u>o| < t- 
The situation just described is abbreviated by the notation 


lim f(z) — u’o. 
o 


Informally, the definition means that (1) the domain contains a disc with positive radius 
and centered at zq , and (2) we may arrange that /(z) lie as close as we want to U'o. 
specifically, within any small distance f, provided that we impose upon z the condition 
that it lie sufficiently close to zq. specifically, within some distance d. The concept of 
limit of complex functions yields the concept of continuity. 


Definition 12. 

but only if. 


A function / : D —* C is continuous at a point z 0 € D C C if. 
lim /(*) = /(*„). 


The function / is continuous on a subset 5 C D if. but only if, the function / is 
continuous at each point r 0 € 5. 


EXAMPLE 16. Let c denote any complex number and consider the function 


C - C. /.(*) = + c. 


Such complex functions will give rise to fractal Julia sets in the next sections. The function 
/ c is continuous at each z 0 € C. as the following argument shows. Firstly, the difference 
/«(*)-/c(*o) factors. 


(1) /((*")" /(( ; o) = (‘* 5 + c) - (:§ t c) = ; 3 - - (: - :o) (i + -o ) ■ 

Secondly, restrict the argument to those points z € C at distance less than one from co¬ 
in other words, assume that |z — zq| < 1. For such points z, the triangle inequality gives 


|r| s= jz — z 0 + ’ol 5k- Zq| + kol < 1 + kol- 
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and. consequently, 

(2) 1 2 + zo| £ jr| + \za\ < (1 + |£ol) + |*oI = 1 + 2|ro|. 
Substituting this result into the factorization (1) of f{z) — /(zq) yields 

(3) i fix) - /(r 0 )l = I*' - *o| • I* + I < \z - zo| • (1 + 2|2 o!) . 


Next, for each positive tolerance (real number) t > 0 let 

d = min 11, -—-—■ 1 
l l+2|£olj 

(the smaller of 1 and f/( 1 + 2|zo|) ). If Jr — ro| < d then \z — rol < d < 1 and the 
estimate (2) holds, and. consequently, (3) holds. Therefore. 

[/(*) — /( ? o)l = I* — *o| • 1* + -ol < !* — *o| ■ {1 + 2Jzo 1) < d{1 4- 2jz 0 |) 


which means that f c is continuous at zo- □ 

The concepts of convergence, limit, and continuity may also be expressed in terms of 
“open sets." as demonstrated below. 

Definition 13. The open disc with radius r and center at z is the set of all 
complex numbers at distance less than r from z, a set denoted by D(r.r): thus. 


D{z,r) = {tr : tr € C and |u* — zj < r}. 

The unit disc is the set £>(0.1), which consists of all complex numbers at distance less 
than one from the origin. 

Besides open discs, w'hich do not include the circle surrounding them, "closed discs." 
which include the circle surrounding them, also arise in the context of Julia sets. 

Definition 14. The closed disc with radius r and cent er z i s the set of all 
complex numbers at distance at most r from r, and denoted by D(r,r): thus. 

D(r,r) = {tr : ti’£ C and |tr — zj < r} . 


Definition 15. A subset U C C is open if, but only if. for each point : f I" there 
exists a positive real number r such that £)(r,r) C V. 

Thus, a subset is open if, but only if. it contains an open disc around each of its points. 

EXAMPLE 17. Every open disc .D(z.r) is an open set. To verify this assertion, 
let w represent any point in D(z,r)\ thus. |u* - r| < r. Next, let s = (r - |tr - z|)/2. 
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Figure 4. {a) Open discs are open sets, (b) The boundary of a 

disc is the circle surrounding it, 

which represents half the distance from tr to the circle surrounding the disc D[z.r). If 
a complex number q 6 C belongs to the open disc D(u'.s). then \q — u-| < s. and. by 
the triangle inequality, 

19- *1 = l5-*r + tr— z\ < to-«’l + K’-i| < (r- |tr- z|)/2 + |u'-r| = r - |ir - zj/2 < r. 

Hence, every point q in £>(u\s) also lies in the disc D{z, r), which means that D[ic.s) C 
D(z.r). as in figure 4a. Therefore, the open disc D(z,r) is an open set. □ 

The concepts of “open disc" and ‘'open set" form the basis for the concept of “bound- 

ary. 

Definition 16. The boundary of a set 5 C C is denoted by dS and consists of 
all points - € C such that for each positive radius r the open disc D(r.r) contains at 
least one point in 5 and at least one point not in 5, 

Informally stated, the definition means that a point ir lies on the boundary of a set 
S if. but only if. there exist points in 5 and points outside S as close as we want to tr. 

EXAMPLE 18. Let 5 = D{0,1) represent the closed unit disc: its boundary dD(0. 1) 
consists of the unit circle, which surrounds it and has radius 1 and center at the origin. 
0. as in figure 4b. 

To verify this assertion, consider any point u' on that circle, at distance 1 from 
the origin: thus, |ir| = |tr — 0| = 1. For each positive radius r consider the open disc 
D(u\r) with radius r and center at the same point u\ Then D(tr.r) contains a point 
p in 5 and a point q not in S, because we may find such points p and q along the 
half-line from the origin through u\ for instance,' p = (1 — [r/2])u\ for which |p| < 1. 
and q = {1 + [r/2])u\ for which |g| > 1. as shown in figure 4b 

Co nsequently, for each posi tive rad ius r, the disc Z>(u'.r) contains a point p in 
D( O.l) and a point q not in D(0,1). Therefore, u' lies on the boundary of the unit 
disc. Conversely, every point on the boundary lies on the unit circle. □ 










24 


Chapter 1 QUADRATIC JULIA SETS 


The following proposition describes the relationship between open sets and continuity, 

Proposition 3. A function f : D —* C is continuous at a point z £ D C C if 
but only if for each open set V containing f(z) there ezuti an open jet U such thot 
z £ V C D and f(U) C V. In particular, if the domain D u itself open, then the 
/unction / u continuous on D if , feut only if. for each open set V C C the inverse 
tmagc / -1 (V*) = {z : z £ D and f(z) £ V} is open. 

Proof. Suppose that for each open set V containing f(z ) there exists an open set 
U such that z £ U C D and such that f(V) C V. For each positive real number t > 0 . 
consider the open set V = D(f(z),t) {the open disc with radius t and center at /(:)). 
Since the corresponding set U described in the hypothesis is open, there exists a positive 
real number c > 0 such that U contains a disc Z?(r,c). Consequently, again by the 
hypothesis, f(D(z,c)) C /(£/) C V = D{f(z).t). which means that if [ir — r| < c then 
w £ D and |/(u') - /(r)| < t. Thus, / is continuous at z. 

Conversely, suppose that / is continuous at a point z £ D. The definition of 
continuity then asserts that for each open set V containing /{r}, which contains an 
open disc of the form D(f(z),t). there exists a positive real number c > 0 such that if 
|u? — z\ < c then u> £ D\ consequently, D{z.c) C D. and |/(u ! ) — /(r)j < t. and hence 
f(D[z.c)) C D{f{z).t) C V. Therefore, the set U — D(z.c) satisfies the requirements of 
the propositions. □ 

The following proposition describes the relationship between continuous functions and 
converging sequences. 

Proposition 4. A /unction f : D —* C w continuous at a point z £ D C C if. but 
only if for each sequence (r*) in C converging to z. the sequence (/(r*)} converges 
to f{z). 

Proof Suppose that / is continuous at z £ D. and consider a sequence (;*) 
converging to z. To show that (f(zt)) converges to f{z). let f>0 represent a positive 
real number. By continuity of / at z, there exists a positive real number c > 0 such 
that D(z,c) C D and f[D{z.c)) C D(f(z),t). By convergence of the sequence (c*) 
to z , there exists an index A' € N such that if k > K then - z\ < c. which 
means that r* £ D(z.c); consequently, /(;*) £ f{D{z.c)) C D{f(z).t). which means 
that |/( 2 t) - f{z )| < t. Therefore, the sequence /(r*) converges to f{z). 

For the converse, proceed by contraposition: assume that / is not continuous at 
z £ D. Then there exists a positive number t > 0 such that for each positive number 
c > 0 there exists a complex number z c with \z c — z\ < c but \f{z c ) - /(;)! > t or 
z c $ D. By induction on k £ N and with c t = 1 /(k + 1), the sequence (r Ct ) converges 
to 2 but f{z Ck ) does not converge to f(z) {some of the terms f{z Ct ) may also fail to 
exist). □ 

REMARK 9. (This remark serves only to explain the title of the present section.) 
The collection T of all open subsets of the complex plane. 


7 = {V : U C C and V is open) . 
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satisfies the following properties, which make T a topology on C: 

(1) The empty set, 0, and the entire plane, C, both belong to the collection: 0 £ 7 
and C € T. 

(2) For each family (set) T C T of open sets, the union of all the open sets in the family 
is again an open set: [JT € T, 

(3) For each finite family Q C T of open sets, the intersection of all the (finitely many) 
open sets in the family is again an open set: f]Q e X. 

3.2, Theorems from complex analysis 


Most concepts familiar from calculus have their counterparts in complex functions, 
for instance, the derivative. 

Definition 17. The derivative at a point z £ D of a function / : D C C —• C 

is the following limit (if it exists), denoted by /'(r): 


/'(.-) = Urn 

h — 0 


f( Z + h)-f(z) 
h 


The function / is called complex differentiable if. but only if, its derivative f'(z) 
exist5 at every point z € D. 

Notice that the definition of the derivative of a complex function coincides formally 
with that of the ordinary derivative of a read function. The complex derivative differs from 
the real derivative only by the fact that all values in the limit ( r, /(r), h. and f(z 4- h) ) 
may be complex; in particular, this means that the complex increment h may tend to 
zero in any fashion in the complex plane Due to the similarities between complex and real 
algebra, and between the topologies on the complex plane and the real line, most formulae 
from ordinary calculus also hold in complex calculus. 

EXAMPLE 19. Consider the function f c : C — C with f c {~) = z 2 4 e. Then 


fM 


lim 
A—o 


f e (z 4 h) - f e (z I 
h 


A—0 h 


= lim 

a— o 


r 2 + 2 hz + h 2 -z 2 
h 


lim (2 z 4 h) = 2 z. 


We shall take the following particular cases of theorems from complex analysis for 
granted. Their proofs require at least a short course in complex analysis, but their state¬ 
ments may be easily understood and are indispensable to the study of Julia sets. 


Theorem 2. (Cauchy’s Ineq uality) Consider a complex differentiable function 
f : D C C —* C and a closed disc D(u\r) C D contained in the domain. Also, let 
M — max{]/(r)| : \z — tr| = r} denote the manmiim value of |/(z)| over all points : 
on the circle with center u' and radium r. Then 
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Theorem 3. (Open Mapping Theorem) Let / : C —* C represent a non¬ 
constant function defined and complex differentiable on the whole plane. Then for each 
open subset U C C the image /(£/) is also open. 

Routine exercise 

Exercise 17. Prove the following assertion 

At a fixed point z m , the following equality holds: (/ ( . on ) , (-.) = (/'(r,))". 


4. PROPERTIES OF QUADRATIC JULIA SETS 

Complex analysis and topology provide indispensable tools to describe the nature of 
"Julia sets." for not only the properties but also the definition of Julia sets rely on topology 

4.1. Two examples of Julia sets 

The topological concepts reviewed in the preceding section provide the language nec¬ 
essary to define the nature of Julia sets, which are the boundaries of filled Julia sets. 

Definition 18. The Julia set of a quadratic polynomial f e with f e (~) — - 2 •+• c 
is the boundary dl\ c of the filled Julia set K c , and it is usually denoted by J c 

EXAMPLE 20. Recall, from example 12. on page 12. that the filled Julia set A’o of 
the squaring function /o ( fa{z) = r 2 ) is the closed unit disc. 

A'o = D(0,1) = {r : :gC and \z\ < 1} . 

Recall also, from example 18, on page 23, that the boundary of the closed unit disc is the 
unit circle; therefore, the Julia set of the squaring function is the unit circle: 

J 0 = <3A'o = {tr : irfC and \w\ = 1} . 


EXAMPLE 21. Consider the quadratic polynomial /_2 : C —* C. with /_ 2 (- ) = 
r 2 — 2. As conjectured in exercise 14, the filled Julia set A’_ 2 of /_ 2 is the segment 
[-2,2], which coincides with its own boundary. Thus. J_ 2 = [—2.2]. To prove this 
assertion, use a composition with the function 


5 : C \ [-2,2] 


C, ?(*) 


0 


(The notation C \ [—2,2] stands for the complement of the interval [-2.2]; thus, the set 
C \ [—2,2] looks like a plane with a slit along [-2.2]). As explained below, the function 
g maps C \ [—2.2] onto the complement of the closed unit disc A’o, in other w'ords. onto 
(r : z £ C and \z\ > l}. as shown in figure 5. The function g has an inverse function. 

: {ir : w 6 C and |tr| > 1} —* C \ [—2, 2], = tr 4- —. 
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-2 



a i—C> g (z) 


I Im 



Figure 5. The function g maps the region outside the interval 
[- 2 , 2 ] onto the region outside the dosed unit disc. 


(The formula tr + 1/u - also extends to the complex numbers it for which |u-] = 1. and 
it explains why g maps C \ [—2,2] onto the complement of the unit disc. If ir = (x. y) 
lies on the unit circle, then |u'| = 1. and 1/u? = W, because u'lr = (x. y )(x. -y ) = 
i 2 + y 2 = |u'j 2 = 1; moreover, [x| = y/x 2 < \fx 2 ■+■ y 2 = |u'| = 1. and. consequently. 
|u’ + l/u'| = |u’ + W\ — |(x.j/) + (x, — y)| = 2|x| < 2. Therefore, if w lies on the unit circle, 
then |u’ 4- l/u'| lies on the interval [—2,2],) 

Straightforward algebra shows that g~*of 0 og = /_ 3 , which means that g ~ 1 (/ 0 (j(c))) 
=/- 2 (r) for every r€C\[-2,2]. Consequently, /° 3 = (ff -1 o f 0 o g)° n s o / 0 en c g 
{because intermediate compositions g~ ] o g cancel out). Therefore, (r„) = (/I 3 (*o)) — 
(ff -1 (/o n (fff*o)))) diverges if. but only if. (u-„) = {f™{g{z 0 ))) does, which means that 
the complement of the filled Julia set of the squaring function. C \ A'o. is the image under 
g of that of /_ 3 . Hence, C \ A'o = g(C \ K- 2 ) and C \ A '_2 = ? _1 (C \ A'o) = C \ [— 2 . 2 j. 
Finally. = <9A '_3 = d[—2,2] = [—2.2] in C. 

4.2. Invariance of Julia sets under their quadratic polynomials 


Except for the particular values c = 0 and c = -2. examined in the preceding 
subsection, the investigation of the Julia set of f c for values of c other than 0 or -2 
requires a more substantia] analysis for each value of c separately. Therefore, instead of 
investigating any particular Julia set, we shall focus on a few general properties of Julia 
sets that rely only upon the analysis and topology reviewed in the preceding section. The 
first property of quadratic Julia sets examined here concerns their invariance under their 
associated polynomial f e . That is, each function f c maps the associated filled Julia set 
A’c onto itself. 

Proposition 5. For each complex constant c € C. the filled Julia set K c is invariant 
under the /unction f e ; that is, f c (K c ) = A' c . 

Proof. For each initial point zq € A' c in the filled Julia set, the sequence of iterations 


"O’ *1 — /c(*o)’ "2 /c{"1 ) * *3 — /c (~ 2 ) ’ 
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remains bounded. Consequently, so does the sequence starting at z\. 

^ 3 i *2 ~ fcM, 2 3 = f cl * 2 )* ^4 ~ fc ( * 3 L .... 

Therefore, z^ = /(z 0 ) also lies in the filled Julia set A' f , which proves that /cfA’d C A\ 
For the converse inclusion, select either preimage of zq, for instance, r_j = ^/z 0 ~ c. 
so that then the sequence that starts from z_] also remains bounded: 


* — 11 ^0 Fe(* — 3 h ^3 — /c(*o)i *2 — f i }, .... 


Consequently, the preimage z_ a also belongs to the filled Julia set A ' C1 which shows that 
A'e C /c(A'c) (because z 0 = / c (*-i) )• □ 

Not only does each filled Julia set A' c remain invariant under the function / c . but 
so does its boundary dK c - J e (the Julia set itself), and its complement, C \ I\ c . which 
has a special name that reflects its properties. 

Definition 19. The basin of attraction of infinity of a quadratic polynomial f c 
is the set, denoted by A^.c, which consists of all complex numbers zq € C that give rise 
to an unbounded sequence. Thus, Aqc.c = C \ A' c - 

The term “basin of attraction of infinity" arises from the fact that, for the particular 
sequences (z n ) -• (/“"(zo)) considered here, if a point zq does not belong to the filled 
Julia set, then the unbounded sequence that it initiates diverges to infinity (recall exercise 
13). Informally, infinity attracts all sequences that start at any point z 0 € . 4 ^ f . The 
following sequence of theoretical exercises shows that each of the sets . 4 ^ c . K c . and J c 
is invariant under mappings by f c . 


Theoretical exercises 

Exercise 18. Prove the following assertion. 

The basin of attraction of infinity. A KiC , u an open set. 

Exercise 19. Prove the following assertion. 

The i>a.«tn of attraction of infinity is invariant under f c . In other words. f c maps 
A x . t onto itself; thus, f e (A x , c ) = A®. c - f~ l {A x , e ). 

Exercise 20 . Demonstrate that / e (J e ) C J c and that f~ 1 {J c ) C J c - 

Exercise 21 . Denote by A‘“ =■ K c \ dh\ the interior of the filled Julia set. 

which consists of all points z € K c that do not lie on its boundary (all points z € K z 

for which there exists a positive radius r > 0 such that D(z.r) c A' c ). Prove that 
/ f (A' c °) = A7 = / c - 1 (A7). 

Exercise 22. Apply the preceding exercises to prove the following assertion. 

For each complex constant c € C the Julia set J c = dK c *3 invariant under the 

associated quadratic polynomial f c and its inverse; in mathematical formulae. f c (Jc) = 
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The property of quadratic Juba sets examined here concerns the presence on the 
boundary J c of at least one of the two fixed points of / c , where the derivative 
has a modulus of at least one. From such a fixed point, the Inverse Iteration Method 
generates only points on the Julia set J c . 

Definition 20. A non-attracting fixed point of a function / : C —* C is a fixed 
point z. where j/'(2.}| > 1. A non-attracting fixed point 2. is called repelling if 
|/'(z.)| > 1 and parabolic (or neutral) if |/'(z.)| - 1. In contrast, at an attracting 
fixed point, |/'{z.)j < 1. (A fixed point where /'(z.) = 0 is also called superattracting. 
but we shall not need so fine a classification.) 

The terminology for the classification of fixed points arises from the behavior of itera¬ 
tions beginning with an initial point near a fixed point. If z 0 lies near an attracting fixed 
point 2,, then the sequence of iterations ( z n ) = (/° n (2o)) converges to 2.. In contrast, 
if the initial point zq lies near a repelling fixed point 2., then the sequence (2 n ) may 
(but need not) diverge away from 2,. 

EXAMPLE 22. With c = 0 the squaring function /o- with fo(z) = z 2 , has two 
fixed points. 0 and 1, since / o (0) = 0 2 = 0 and /o(l) = l 2 = 1. The first fixed point is 
attracting, whereas the second fixed point is repelling, as the following arguments confirm 

(22.1) At 2 = 0, / 0 (0) = 0 and |/'(0)| = |2 x 0[ = 0 < 1. 

(22.2) At 2 = 1, /o(l) = 1 and |/'(l)j = |2 x 1| = 2 > 1. 

Moreover, from examples 12 and 20, recall that the filled Julia set A’o consists of the 
closed unit disc. A'o = D(0,1). and that its boundary, dKa = Jo, is the unit circle 
Jo = {r : z € C and \z\ — 1}. The attracting fixed point 0 £ A'o belongs to the filled 
Julia set, whereas the repelling fixed point 1 £ J 0 lies on the Julia set. 

EXAMPLE 23. With c = 7* the function fi/ A : C —* C. /i/ 4 (*) = z 2 4- V* has 

only one fixed point, at 2, = V21 shown in example 13. on page 16. This fixed point is 

parabohc, because |/X<( 72)1 = |2 x 72! = 1- This particular parabolic fixed point lies on 
the Julia set J-y / 4 = 3A*i/ 4 . as the following proposition shows. 

Proposition 6, The parabolic fixed point 2. = 72 he$ on the Julia set J]/ 4 . 

Proof. Recall that proposition 2, page 17, (with n = 0 ) established that every fixed 
point lies in the filled Julia set; consequently. 72 € A‘i/ 4 . To show that 72 lies on the 
boundary, it suffices to prove that each open disc D( 72. r) contains a point z 0 not in 
A"i/ 4 . To this end, consider the point z 0 = 72 + jr/2], which initiates the sequence 

z o = 72 + [r/ 2 ], 

2i = z 2 + c = ( 7 2 + [r/2]) 2 + 7, = 7 2 + (r/2)(l + r/2)„ 


Comparing the distances from 72 to *o and 23. observe that 
* 1 - 2 . 72 + (r/2)(l + r/2)- 72 (r/2j(l 4- r/2) 


*0 — - • 


72 + [r/2] - 72 


r/2 


= 1 + r/2 > 1. 






30 


Chapter 1 QUADRATIC JULIA SETS 


Hence, induction shows that |z n — z 0 | > (1 + T l 2)”, which diverges to infinity. C 

In contrast to the situation with c — */*, for every c ^ J /» the quadratic polynomial 
fc has at least one repelling fixed point, and even, 1 such fixed point belongs to the boundary, 
as the following tw-o propositions show. 

Proposition 7. If c ^ '/< then f c has a repelling fixed point. 

Proof. Recall from example 14 that for c ^ 7< the quadratic polynomial / r has 
tw'o distinct fixed points, expressed by the quadratic formula 

1 ± \/l ~ 4c 
2± = 2 ' 

At the fixed points, _ 

/e(*±) = 2r± = 1 ± 7l - 4c. 

If either fixed point, for instance r + = (1 + y/T^~£c)/2< is not repelling, then 1/1(2+)| = 
jl 4- \fl - 4c| < 1. For the rest of the argument, recall from definition 1 and remark 3 that 

M = |(*,y)l = V* 2 + V 2 ~ \/Re(;) 2 + Im(z) 2 , 

In particular, if y — Imfz) ^ 0 then |z( > |ij. Thus, Re{ v^l — 4c) < 0: since c ^ '/<• it 
follows that l-4c ^ 0. and if Ref \/\ - 4c) were zero then Imf^/l - 4c) would not be zero, 
and the Pythagorean theorem would show that |1 + y/l — 4cj > 1. contrary to the initial 
assumptions. With Ref — 4c) < 0 it follow's that Re(l 4 — 4c) < 1. which means 

that 1 + \/l — 4c lies in the left-hand half-plane bounded on its right by the vertical line 
through x — 1. Consequently, at the other fixed point. Ref —\/l — 4c) = -Ref \/l - 4c) > 
0 and Ref] - y/l - 4cJ > 1. Therefore, |1 - \/l - 4cj > jRefl - \/l - 4c)| > 1. which 
means that the other fixed point, z_, is repelling. □ 

Proposition 8. Every repelling fixed point belongs to the Julia set. 

Proof. By contraposition, suppose that a fixed point z. does not belong to the Julia 
set J c = dK c . Since every' fixed point lies in the filled Julia set A' c . it follows that there 
exists an open disc D(z, ,r) C A' c . (If there were no such disc. then, by definition of the 
boundary, z. would lie on the boundary J c .) Also, since the complement A XiC — C\K c 
is an open set, it follows that Z?{z.,r) C A' c because no point on the boundary of the 
disc D(z.,r) may lie in A 0CiC . Since the filled Julia set remains invariant under the 
function f c , it follows that / c on (A' c ) C A' e C D(0, A/), and hence it follows from Cauchy's 
inequality and exercise 1? that 

)r = K/r )’(--.)i <^. 

Thus, the powers {f' c (:,)) n remain bounded, w'hich has the necessary consequence that 
^ 1- Therefore, 2, is not a repelling fixed point. □ 

Corollary 1. All the preimages of every repelling fixed point of f c belong to the 
Julia set J c . 
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Proof. Every repelling fixed point of f c belongs to the Julia set J c . and all preimages 
of all points on the Julia set remain on the Julia set { / C -] (J C ) = J c ) by exercise 22. □ 

Though all the preimages of a repelling fixed point belong to the Julia set, there might 
not exist enough of them to outline the Julia set. Fortunately, the following proposition 
guarantees that even,’ non-attracting fixed point has infinitely many distinct preimages. 

Proposition 9. For each complex constant c € C every non-attracting fixed point 
2, °f fc has infinitely many distinct pretmages. 

Proof. Observe that if the preimage of a fixed point 2 . consists of 

only one point, namely / c _1 ({2.}} = { 2 .}, then z m is a double root of the equation 
z 2 + c ~ x. , which has the necessary consequence that f^(z.) = 0, and. consequently, 
2 z. — 0. Therefore, 2 . is an attracting fixed point, which contradicts the assumption. 

Thus, a non-attracting fixed point has at least two distinct preimages, which we denote 
by 2 , itself and z L Also, denote by Z the set of all preimages of 2 .; thus. 

2 = 0 (A")*’«*.})• 

nss J 

To obtain a contradiction, assume that Z contains only finitely many preimages, and 
let A" represent the smallest number of backward iterations necessary to produce all the 
preimages: 

N = min < m : m 6 N and Z = ({*.}) 

L n = l 

Notice that „V > 2 because / c -1 ({*.}) ~ {r.. 2 t } contains a preimage different 
from the fixed point r,, and because /'*({;*}) contains neither 2 , nor r f , since 
/ c (:.) = z. / z i and f e ( r f ) = 2 . ^ Next, let 



N -1 

e U (/D" ({*•)) 

ns 3 

represent a point obtained after exactly N back w-ard iterations of f c . Thus, f° s ( 20 ) = 2 . 
but f° n (zo) ^ 2 , for every n < A\ Moreover, select any preimage of 20 , and denote it 
by r_] € /“ 1 {{* 0 })- Then 2 _] is itself a preimage of 2 .. because 

/U +1 (.--i > = r, K 

Consequently, there exists an integer m € {2. A"} such that 2 _] € (/° m ) _1 <{-.}*' 

and, therefore, 

2 . = /;-(*_,)■ /;— 1 (f c (:-i )) = /c em_1 (‘o), 

which has the necessary consequence that z 0 € * ({;.}), contradicting the hy¬ 

pothesis on zq. □ 

The foregoing results guarantee that if the Fixed-Point Inverse Iteration Method starts 
from a repelling fixed point, or from the parabolic fixed point if c = 1 /a. and if it retains 
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all the preimages that it computes, then it generates infinitely many points, all of which 
belong to the Julia set. In practice, however, many implementations retain only one 
preimage 2 _„ at each stage and then select the next preimage z_( n+ u = ± J~z _„ - c 
through a random choice of the sign ±. as in the program in exhibit 4. Yet such a random 
choice does not ensure {but only malces it improbable) that the algorithm will not remain 
trapped in a loop through finitely many preimages. A method for selecting the sign ± at 
each stage so that the algorithm generates a sequence of infinitely many distinct preimages 
constitutes the topic of the second research problem in this chapter. 


Non-Attracting Fixed-Point Inverse Iteration Method. The fol¬ 
lowing modification of the Fixed-Point Inverse Iteration Method generates only 
points on the Julia set. 

Step V, Compute 1 + — 4c. If |1 + \/l — 4c| > 1, select the fixed point 


*• = *+ = 


1 + y/l ~ 4c 


Otherwise select the other fixed point, which is 

1 - \/l - 4c 




Next, proceed exactly as in the Fixed-Point Inverse Iteration Method, selecting 
either square root randomly at each step. 


REMARK 10. There exists a theorem stronger than proposition 9 (which guaranteed 
that non-attracting fixed points have infinitely many preimages). Specifically, the set Z 
of all preimages of a non-attracting fixed point r, on the Julia set, as described in 
proposition 9. is dense in the Julia set. Informally, this means that there exist preimages 
r_ n of 2 . as close as we want to any point z € J c on the Julia set. Mathematically, this 
means that for each point z € J c on the Julia set. and for each open disc Z?(z,r) about 
2 . there exists a preimage ?_„ €in D{z, r) H J c . For a proof of this assertion, 
which relies upon Montel’s theorem about normal families of complex analytic functions, 
consult Linda Keen's chapter, “Julia Sets," 1 page 62. in Devanev and Keen's book [32]. 

REMARK 11. Though the Non-Attracting Fixed-Point Inverse Iteration Method 
may generate an infinite and dense set of points on Julia sets, comparisons with pictures 
obtained through different algorithms reveal that such points are not evenly distributed 
on the Julia set. As a consequence, the resulting plots show many points in some parts 
of the Julia sets, and very few in other parts, which yields a poor resolution of the finer 
details of Julia sets. For other algorithms, consult the books by Peitgen et al. [36]. [37]. 
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Julia sets 
on the HP-28C&S 


« c RBS 

4 * 

1 *- T 1 

2 

^ ' r- • STO r- 

32 ’ H 

* STO 

H 137- * 

2 x 

* s * STO 


* r- R -*■ C 

PMRX 


Sii^ g » 0 eKt£r 


• *e t UP ' STO 


« 1 c 4 

* — 

T 1 < 1,0 > 

* DUP 


I F RBS 

1 < 

THEN 2 

NEG 

— 


e8V£r 

' f i xeO 

■ STO 

* CLLCD 

1 2000 

STRRT c 

- T 

< 1 • 0 > 

RRND 

IF .5 

< 

THEN 

NEG 

END * 

PIXEL 

DUP 

£?§g l £ R 

£Rter 

' beck 

* STO 


« setup 
f i * e c) back * 
ti NEG R-*C 
PMRX c NEG r- 
t-i •+• NEG R -► C 

ESifc 1 t^SfER 

* Jul i a • STO 

.25 * c 1 STO 

USER Julia 

n 

u 


Comments on the programs 
for Julia sets 


Julia sets 
on the HP-48SX 


Subroutine setup computes a 
bound for the filled Julia set, 

1 + JT+ 4|c| 
r = r, — ---, 

C 2 

and adjusts the screen so that it 
covers the square [—r, r] x [-r, r], 
hence also the filled Julia set; 
h equals the width of each pixel. 
ENTER and STOre in 'setup' 

Subroutine fixed calculates anon- 
attracting fixed point: for the nu¬ 
merator. if |1 + i/I - 4c! < 1 
t hen 2 - (1 + y/1 - 4c = 1 - 
\/l — 4c and |1 - \/l - 4c| > 1 
if c ^ explained in the 

text). 

ENTER and STOre in 'fixed' 


For the HP-28C&:5 only, subrou¬ 
tine back iterates the computa¬ 
tion of preimages, 

/-■({.-}) 3 ±V7=~c, 

with the commands c — and 
with a random choice of sign. 
ENTER and STDre in 'back' 


The main program. Julia, com¬ 
putes and plots preimages of } e ( z ) 
= z 1 + c {alone on the HP-48SX, 
with back on the HP-28C&:S). 

ENTER and STOre in 'Julia' 

Tutorial: 

STDre a complex number in 'c' 
execute Julia 
in the USER or VAR menu. 


« a BBS 4 * 

1 ♦ 4" 1 ♦ 2 

/ ' r- * STO r- 

64 ^ 1 h> 1 STO 

ERASE r 128 

x 13 1 * r- 

R -*• C DUP NEG 

SWAP PDIM # 

1 3 1d # 1280 

PD IM »■ ENTER 
• set up 1 STO 


« 1 C 4 * - 

T 1 < 1,0 > 

* DUP 

IF RBS 1 ^ 

THEN 2 - 

NEG 

END 2 / DUP 

PIXON» ENTER 
* f i xeO 1 STO 


« s et up 
f i x ed 1 2000 

STRRT C-T 
< 1 • 0 > RRND 

IF .5 < 

THEN NEG 

END * DUP 

DUP NEG 
PIXON PIXON 
NEXT DROP 
* ENTER 

'Julia' STO 

. 25 * c ' STO 

VRR Julia 



Exhibit 4. 
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Exercises 


Routine exercises 

Exercise 23. Verify that /_ 3 y 4 : C —♦ C. /_ 3 y 4 (z) ~ : 2 ~ 3 /4 has a parabolic fixed 
point at z_ = —’/ 2 and a repelling fixed point at — 3 / 2 . 

Exercise 24. Identify the nature (repelling, parabolic, or attracting) of each fixed 
point of the function /_ 2 : C — C. with /_ 2 (r) = z 2 — 2. 

Theoretical exercises 

Exercise 25. As the preceding examples and exercises reveal, not all quadratic 
polynomials of the type / c (r) = z 2 ■+ c have a parabolic fixed point. In fact. j e has 
a parabolic fixed point if. but only if, the constant c lies on a particular curve in the 
complex plane. Identify that curve. {Write the two equations that mean that f c has a 
parabolic fixed point. Then transform the resulting equations into one polar equation, and 
identify the polar equation with that of a standard curve.) 

Exercise 26. Prove that an attracting fixed point z. of f c cannot lie on the Julia 
set. but belongs to the interior A'° of the filled Julia set. Moreover, if z 0 lies close 
enough to then (z n ) tends to z % . (The proof may require a complex line integral.) 

Experiments and conjectures 

Exercise 27. Try the Non-Attracting Fixed-Point Inverse Iteration Method and 
generate pictures of the Julia sets J c corresponding to the following constants: 

(27.1) c = (-0.12256117,0,74586177) This constant c satisfies the equation 

c 3 + 2c 2 + c + 1 = 0 (which yields theoretical properties of the Julia set 

J c ) and J c is called “Douady's rabbit;” see Devaney's An Introduction to 

Chaotic Dynamical Systems, 2nd ed., [30], plate 3. 

(27.2) c = (0.360284, 0.100376) J t is called “the dragon;’* see loc. cit. plate 1. 

(27.3) c = -1.54368901269 This constant c satisfies the equation c 4 1 - 4c 3 + 

6c 2 + 6c+ 4 = 0. and J c is Devaney’s dendrite: see loc. ext., pages 293-295. 

(27.4) c = t = (0,1) The Julia set J x is called “the dendrite,” 

(27.5) c = 2 /s The resulting Julia set J 2 / 5 consists of many disconnected pieces. 

The following two exercises investigate the results of the Fixed-Point Inverse Iteration 
Method with a predetermined sequence of signs to select either square root at each step. 

Exercise 28. Considering the squaring function, with / 0 (z) = z 2 , start from the 
fixed point = 1, select z_j = — = — 1, and. henceforth, always choose the square 

root described in exercise 5 (which coincides with the square-root key on supercalculators): 

r 0 = z + = 1, r_! = —1, z-i = = i = (0,1), 2-3 * V / tM) ■ (1/V5, i/VS). ... 

Compute and plot several additional preimages, and describe the result. 
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Exercise 29. With the function f l / 4 : C — C, fj/ 4 (z) = z 2 + >/,. start from it? 
unique fixed point, 2 . = 1 / 2 . and alternate the choice of the sign ( - or + ). which yields 

‘0 — = l h, 

~-l = ~V Z 0 - c = - yj l h — 'A = — "s/^A = — V 3 * 

* — 2 = c ” + V 2 ~ = +\/— 3/4 = iv^/2. 

*-3 - - c= -^/iv^/2- *A= -^/{-V^v^/S) 

— (—0.570695986873 ..., —0.758744956776 ...), 


Perform about two dozen iterations and describe the result. 

Exercise 30. How do Julia sets relate to the Butterfly Effect in chaos 0 Answer 
this question with a short essay based upon the properties of Julia sets developed in the 
present chapter and upon the description of the Butterfly Effect in the introduction. (Of 
course, you may also consult other references, but doing so is not required.) 


5. EXPLORATORY TERM PROJECTS 

The present section proposes two projects that may be attempted in parallel with the 
study of this chapter (hence the phrase “term projects”). Since the results involve only 
elementary complex variables and are known to specialists, the projects do not qualify as 
research. However, some of the details are not readily available in many sources, which 
means that you may have to search the literature for some results and derive others yourself 
(hence the adjective "exploratory"). 

5.1. Julia sets of general quadratic polynomials 

Project 1 . Instead of the particular type of quadratic polynomial considered so far. 
f c wuth / t (:)r: ! + c, consider every quadratic polynomial /. for which 

/(«’) = pu ’ 2 + qw+j, 

with complex coefficients p, 9 , and j, such that p ^ 0 . 

How do the Julia sets of all quadratic polynomials / relate to the Julia sets of the 
particular type J c for f c ? 

To arrive at a conjecture (a tentative answer), you may either modify the algorithms 
given in this chapter and then plot many Julia sets for various p, q. and or you may 
proceed theoretically as in example 21 , on page 26. In either case, however, you then ought 
to prove your conjecture (or revise it and start again). 
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5.2. A better algorithm for Mandelbrot’s set 

Project 2. Instead of fixing c£C and examining what sequences (r„) emerge from 
each point r 0 . fix -o = 0 and for each c £ C consider the sequence (tr„) — (/ c sn iOj|: 

0. / c {0) = 0 2 + c= c, / C (c) = c 2 +c, f e (c 2 + c} = (c 2 + c) 2 + c. 

The Mandelbrot set is the set, denoted by ,Vt, that consists of all complex numbers 
c £ C for which the sequence (tr„) = (/° n (0)) remains bounded. 

(1) Verify that —2 £ .Vt, 

(2) Prove that if c £ M then |c| < 2. 

(3) Prove that Mandelbrot's set is symmetric with respect to the horizontal axis. 

(4) Modify and test the "Direct Iteration Method" to plot the Mandelbrot set M. in a 
fashion similar to that for filled Julia sets. 

(5) Prove that your algorithm is better than the ones found in the literature. 

6. RESEARCH PROBLEMS 

The two research problems proposed here concern two problems already mentioned 
in the text and investigated in the exercises. As is the case for all mathematical retearch 
problems, the literature does not appear to offer any solution to these problems, and 
nobody knows in advance what methods may lead to solutions. 

6.1. An optimal bound for quadratic Julia sets 

As proved in the text, no point z £ K c in the filled Julia set of the quadratic 
polynomial defined by f c {z) = z 2 4- c may have a modulus greater than R c = max{2. |c|}. 
Vet the bound R c is unnecessarily large for most Julia sets, because, as prov ed in the 
exercises, no point z £ A’ c may have a modulus greater than r c = (1 + yf \ + 4|c|)/2 < R c - 
Nevertheless, the better bound r c may also be unnecessarily large. Find a better estimate 
of the size of the filled Julia set. 

6.2. An infinite sequence of preimages 

Recall that all the infinitely many preimages of every repelling fixed point — or of the 
parabolic fixed point if c = % — of the quadratic polynomial f c lie on the Julia set 
itself (the boundary of the filled Julia set), and that therefore, keeping all the preimages 
at each step of the Fixed-Point Inverse Iteration Method generates infinitely many points 
that outline the Julia sets J c . Yet also recall, from exercises 15 and 16. that plotting only 
one preimage at each step may trap the Fixed-Point Inverse Iteration Method in a loop 
through only very few points scattered on the Julia set without providing any satisfactory 
picture of J c . Find a way to select either square root { 4- or - ) at each step so that the 
Fixed-Point Inverse Iteration Method does not generate any point more than once. Prove 
that your algorithm satisfies the requirement just stated. 
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7.1, Solutions to the exercises 


Exercise 1. 

(M) (0.0) + (2, 3) = (2,3) 

(1.2) (2,3)+(0,0)= (2, 3) 

(1.3) (1.0) x (2.3) = (2.3) 

(1.4) (2,3) x (1.0) = (2,3) 

(1.5) (2.3) + (4. 5) = (6, 8) 

(1.6) (4.5)+ (2,3) = (6,8) 

(1.7) (2.3) x (4.5) = (-7,22) 
(l.S) (4.5) x (2.3) = (-7,22) 
(1.9) (2.3)+ (-2,-3) = (0,0) 

(1.10) (-2,-3)+ (2.3) = (0.0) 

(1.11) (2.3)(2/j 3, — 3 /n) — (1,0) 


(1-12) ( 2 /i3,-Vn)(2,3) = (1.0) 

(1.13) ((2,3) + (4,5)) + (6.7) = (12.15) 

(1.14) (2.3) + ((4. 5) + (6, 7)) = (12.15) 

(1.15) ((2, 3)(4, 5))(6, 7) = ( — 196, S3) 

(1.16) (2, 3)((4. 5)(6. 7)) = (-196.83) 

(1.17) (2,3)((4,5) + (6.7)) = (-16,54) 

(1.18) ((2.3)(4.5)) + ((2,3)(6,7)) = (-16.54) 

(1.19) (-1,0) x (2.3) = (-2.-3) 

(1.20) (2,3) x (-1.0) = (-2.-3) 

(1.21) (1,0) x (1,0) = (1.0) 

(1.22) (0,1) x (0,1) = (-1,0) 


Exercise 2. (2.1) 

x 5 = \z\ 2 = x 2 + y 2 . 


z + ; = 2Re(r) = 2x. (2.2) z - z = 2lm(z)i= 2 yi. (2.3 i 


Exercise 3. |r 2 [ = \zz\ = \z\ ■ |r| - \z\ 2 . Arg(rr) = Arg(r) + Arg(:) = 2Arg(z). 

Exercise 4. |z| = |z — w + u>j < \z — u’| + |twj, hence |z[ — ]u'| < \z — trj. If |r| > [tr| 
then \\z\ — |u'|| = \z\ - |irj < \z — u>|. If \z\ < |u’| then permute the roles of r and u\ 

ExerciseS. Observe that if w = (u, v) then w 2 = (u, t’)(ti. r) = (u 2 — t’ 3 , ur + ru) = 
(u 2 — r 2 ,2ur): thus. 


w 2 = (u.u ) 2 = 



with 


2 2 

= li — V — 



and 


, X + y/x 2 + y 2 . / — x + yfx* + y 2 

2 uv = 2 \l -r-— sign(y )< 
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Hence. 



X + y/x* ■+• y 2 -x + x/x 2 + y 2 



x 2 + y 2 - 


—j 


= (x.y! = c 


Exercise 6. ( 6 . 1 ) (—u.’) 2 = u ’ 2 = z. (6.2) All square roots 9 € C of c satisfy 

the quadratic equation q 2 - z = 0, and such a polynomial equation, of degree two in q. 
cannot have more than two roots. 

Exercise 7. Apply the result of exercise 3: since u? 2 = z, it follows that |r[ - |ir 2 — 
ju-| 2 . hence |u>| s= and Arg(z) = Arg(u' 2 ) = 2Arg(u’), hence Arg (ir) = Arg (: )/2. 

Exercise 8. (8.1) = ±(0,1) = ±1. (8.2) y/^4 = ±(0,2) = ±2i. (S.3i 

v/(0Tj= ±(1/V2A/V2). 

Exercise 9. Complete the square, verbatim as in precalculus: 



0 — az 2 + bz + c = a 


Hence, dividing both sides by a. rearranging terms, and taking square roots yield 



Exercise 10, r + = (-3,4) and *= (-3, -4); z ~ (—6 ± i/6 2 — 4 - 25)/2. 
Exercise 11. If z 0 6 K c then — 20 £ A' e because Z] = z£ + c = ( — z 0 ) 2 -t- c and 


the sequences starting at z 0 and — z 0 coincide from z\ and beyond. 
Exercise 12. Set d = |zo| - r c > 0, so that (zo| ~ r c + d. Then 



Use the expression just found and the reverse triangle inequality to obtain 


MH*J-e|>||* o r-|c||> 



To complete the induction, assume that 
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holds for some integer n £ N‘, and repeat the preceding argument with d replaced by 
d( 1 + v/1 + 4|c|)\ which leads to the corresponding inequality for n + 1 . Therefore, if 
ir 0 |>{i +ymr c |)/2 then (:„) diverges to infinity. 

(12.2) Start the sequence at z* instead of at z 0 . 

(12.3) If |c| < 2 then r e = (1 + y/l + 4[c|)/2 < (1 + y/1 + 4^2)/2 = 2. If c = 2 
then r c = 2 = |c|. If |c| > 2 then jc| — 1 > 1, hence (2|c| — 1 ) 2 = 4 !c[ 2 - 
4jc| + 1 = 4|cj(|c| - 1) + 1 > 4|c| + 1; consequently, r e = (1 + y/l + 4jcj j/2 < 

(1 + y^( 2 |c| - l) ! )/2 = |c|. Therefore. r c < max{ 2 .|c|} for every complex 
constant c. 

(12.4) If c €] “ oc, 0 ] then — c — |c] and 1 — 4c = 1 + 4|cj. Consequently. 
z+ — (1 + 71 — 4e)/2 = (1 + \fl 4- 4|c|)/2 = r c , which shows that one of the 
fixed points, z+, lies at distance r e from the origin, but that each fixed 
point belongs to the filled Julia set, and so does its opposite — z + . by the 
preceding exercise. 

Exercise 13. If (z n ) is unbounded then there exists an index .V £ \ such that 
kv|>r c = (l + v/TT 4|cj)/2, whence (z„) diverges to infinity*. 

Exercise 14. 

(14.1) z + = 2 and z_ = —1: solve /_ 2 (*) = z, z 2 — 2 = z, z 2 — z — 2 = 0, 

* = (l±v/l-4(-2))/2. 

(14.2) From z+ = 2 , z_, = -y/2 -(- 2 ) = - 2 , z _ 2 = y/-2- (- 2 ) = 0 , z _ 3 = 
± v /0-(-2) = ±v / 2- 

From z_ = — 1 , z_j = 1 — (- 2 ) = 1 , z _ 2 = ± y/i — (— 2 ) = ±73. 

r _3 = ±.y/ 75 + 2 . 

(14.3) If -2 < z_ n < 2 then |z. (n+1) | = y/z-(-2) < y/2 - 2) = 2 . 

Exercise 15. 

(15.1) z+ = 1 and z_ = 0: solve /o(z) = z, z 2 = z, r(z — 1) = 0. 

(15.2) From z_ = 0 , the inverse iteration method remains trapped at zero: if 

z_„ s 0 then z_(n 4 3 ) = y/0 + 0 = 0 . From z + = 1 , the inverse 

iteration method generates preimages z_„ that are all distinct and all lie 
on the circle with radius 1 and center at 0 . 

Exercise 16. 

(16.1) z+ = 3 / 2 and *- = - l h- solve /_ 3 / 4 (z) = z, z 2 - 3 / 4 = z, z 2 - z- 3 / 4 = 0 - 

(16.2) z_ = -V 2 , z_i = 1 / 2 , *-2 = 75/2, z _3 = (^3 + 275}/2, z_ 4 = 1.454.... 

z _ 5 = 1.4S4 ..z _ 6 = 1.494 ..., z _7 = 1.498 ..thus, the sequence appears 
to converge to 1.5. 
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(16.3) Since r_( n+] ) = - v /r_ n + 3 / 4 , consider the function 
9 ■ [ J A- 3 /a) —* ['A- 3 /a ], ff(x) = \/x + 3 A- Since 
g'(x) = 1 /( 2 v^T^), it follows that < 

|$'{ >A)| = 72 < 1 for even* x € [ 74 . 3 /x]■ By the real 
Mean Value Theorem, if z n C [74-72] then there 
exists a real number x € [ 7 -t. 3 A] such that 

*-<»+!) ~ 72 = 9(z- n ) - 5(72) = (r—« - 7a) 



Consequently, |*_ (b+1) - V 2 I < \z~n - 721/2. and, 
therefore, the sequence of preimages (r_ n ) converges 
to := 7z' Exhibit 5. 

(16.4) From z+ = 72. the sequence of preimages (x_„) traces exhibit 5. 

Exercise 17. Proof. Apply the chain rule. (/ 03 )' = (/' 0 g)g': 


U: 2 )' (r.) = ft (Ms.))ftM = /'(*.)£(*.) = U' c (z,)) 2 . 


Then apply induction on n. □ 

Exercise 18. Proof. The results of exercise 12 establish that a point initiates a 
divergent sequence if. but only if, there exists an index n € N such that \:„\ > r c . T his 
means that some term z n lies outside the closed disc with radius r c = (1 + y/T+W\)/% 
thus. r n € C \ D(0.r c ). Since z n = f° n (z 0 ), this also means that 

^ = (/rr 1 (*.) s (rr 1 (c \wyri). 


Consequent ly, 

oc 

Aoc.c = \j (frr 1 (c \D{ o.r c )) , 

n = 0 

Since C \ D(0,r c ) is am open set, and since /*" is continuous, it follows that each 
preimage (/ c on ) -1 \ £)(0,r e )^ isopen. As a union of open sets, A x ^ isopen, □ 

Exercise 19. Proof. To show that /cfAsc.e) C Aoc, c . suppose that :q € .4;x, c . 
Then (z n ) diverges to infinity, which means that the sequence that starts at rj = /( co) 
also diverges to infinity, which also means that z ] € .4oo.c- Therefore, f c {A-x.e) C A x . c . 
and. equivalently, 4,c C / t ''(4,t); thus, f c {A^ rC ) C A«.e C / e _1 (^». e ). 

For the reverse inclusions, select either preimage of ro. for instance r_! = 7*o — c - so 
that /c(r-i) = xo- If z 0 € Aoc.c then the sequence (c„) diverges to infinity, and so does 
the sequence that starts from Z-\ since both sequences have the same terms with indices 
shifted by one. Consequently, € Aao.e, which means that A a 0 , c C /cfAac.c). Since the 
same argument holds for all preimages of zq € it follows that /“'(Aoc.e) C .4^ if . 

□ 
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Exercise 20. If zo € fc(Jc) then there exists a point r_j g J c such that 
/ c (z_i) = zq. By continuity of / c , for each open disc Z?fro.r) the set f~ 1 ( D( ro. r)) 
is open and contains r_i; hence it contains an open disc D{c_i,s) C /“* ( D( r } j. 
Since r_i g J t = dK c , it follows that Z?(z_ t .s) contains a point p g A' e and a 
point q g .4 Xi c. Consequently, f c (p) g f e (K e ) O f c (D{:-i . $)) = A' c nl)/: 0 .r) and 
fc(q) € /cf.Aac.e) n / c (£)(r_i,s)) - nD(: 0 .r). Therefore. z 0 g dh' c . 

Moreover, the preceding argument also shows that if z_i g J c then / c (z_i) g J c . 
which proves that J c C f~ 1 (J e ). 

Exercise 21. By proposition 5, f c (K°) C / c (A' e ) C A' c . and, by the Open Mapping 
Theorem. / C (A7) is an open set. Hence, f c {K°) C A'®. If also follows that A"“ C 
/ c -1 (A’°)- Similarly, / C -1 (A'°) C /ci(A' c ) C A' c by proposition 5, and f~ l ( A'®} is open 
by continuity of / c . hence f ~ 1 ( K °) C A"°. 

Moreover, if z 0 £ A'° then there exist preimages of z 0 - For each such preimage. 
--1 £ .4^ c by exercise 19. and z_j £ 5A' e = J c by exercise 20. Therefore. r_! g 
C\(.4 oc . e UJe)*= A7 and A'° C f e (K°). 

Exercise 22. If zq g J c then there exist two preimages z+ and z_ (which 
need not be distinct from each other) such that / c (z + ) = z 0 = / c (z_). provided by the 
quadratic formula (or the Fundamental Theorem of Algebra). If zj denotes any preimage 
of Zq (thus. z_j g {z+,z_} ), then z_i £ K° by exercise 21, and z_ t £ by 

exercise 19. Therefore. z_i € C \ (A*“ U A^. c ) = 3A' C = J c . and z 0 = / c (z-i) € /c(A c ). 
hence J c C / c (A e ). 

In addition, the preceding argument also shows that if z_] g then r_i € 

which proves that f~ l {J c ) C Jc 

Finally, the results just obtained in the present exercise, and those from exercises 20. 
establish that fdJc) — Jc — / C -1 (A C ). 

Exercise 23. )/l 3/4 (— 1 A) = |2 x (- »/ 2 )| = [ - 1| = 1 and |/l 3/4 ( 3 /*)l = |2 x 3 / 2 | = 
3 > 1. 


Exercise 24. Both are repelling. —1)| = j2 x ( —1)| s 2 > 1 and 1/1 3 (2)| = 
|2 x 2| ■ 4 > 1. 


Exercise 25. 

Solution with conformal mappings. |1 ± \/l — 4cj = 1 
if. but only if, 1 ± \/l — 4c lies on the unit circle, which 
means that ±\/l — 4c lies on the circle C with radius 
1 and center at (1,0). Hence, the squaring function maps 
±s/\ - 4c to 1—4c but it also maps the circle C onto the 
cardioid with polar equation r s 1 + cos{0) (see Spiegel’s 
Shaum’s Outline [14], chapter 8, page 210, mapping C-2). 
A rotation by one-half turn (a multiplication by — 1 ) maps 
1 — 4c to 4c — 1, then a translation to the right by 1 gives 
4c. and finally a compression by a factor of */< yields c 
and the cardioid in exhibit 6. 



Exhibit 6. 
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Solution with algebra in polar coordinates. 1 ± \/l — 4c = (cos 6. sin 8) = t' 6 . Thus. 

1 - 4c = (e ,e - l) 2 = e l2S - 2e ,s + 1 = (cos(20) - 2cos0 + 1 , sin{20) - 2 sin 8, 

s (2(cos 9) 2 — 2 cos 9 , 2 sin B cos 9 — 2 sin 0) = 2(cos 0 — 1 )(cos 9 . sin 9) 

which has modulus r(0) = 2(cos0 — 1) and argument Arg (cos 9, sin 8) = 9. Therefore. 
1 - 4c lies on the cardioid with polar equation r = 2(cos 8 - 1). 

Exercise 26. If is an attracting fixed point of f c then j/'(r, )| < 1. By continuity 
of /’, there exists an open disc D(z m , r) such that |/'(z 0 )| < l fi -I- |/'{r.)j/2 < 1 for 
each 2 q € D( 2 .,r). With an integral along the straight line segment from c, to r 0 . the 
Fundamental Theorem of Complex Calculus yields 

1*1 - *.| - |/«t.-o) - /«(*.)! = \j'° 

< I--0 - --.]sup{|/;(--)| ■ 6 r)} < |.-o - r.| ■ (1 + |/;(.-.)|)/2 < 

Hence, induction shows that (r„) converges to r.; in particular. z 0 6 Ah. 

Exercise 27. 

(27.1) See exhibit la, in the introduction to this chapter. 

(27.2) See exhibit lb. 

(27.3) See exhibit lc. 

(27.4) See exhibit Id, 

(27.5) See exhibit le. 

Exercise 28. The sequence of preimages clusters near the points -1 and 1. 

Exercise 29. The sequence gets trapped and alternates between two points. 

Exercise 30. The Butterfly Effect consists of the great sensitivity of changes in the 
weather due to small initial perturbations; such sensitivity may give us the impression of 
chaos, because we do not know what small initial perturbations cause the large changes 
that we observe. Similarly, sequences (z n ) = (/° n (ro)) of iterations of a function f c may 
exhibit great sensitivity to the initial point, 2 q, especially if zo l' es ° n the Julia set J c : 
every perturbation that keeps zo in the filled Julia set K c keeps the sequence bounded, 
whereas any perturbation of z Q into the complement C \ A' r of the filled Julia set — 
even one too small to be measured — causes the sequence to diverge to infinity. 

7.2. Hints for the term projects 


Project 1. The present line of thoughts expands on a comment in Falconer's Fractal 
Geometry, [34], page 204. and shows that the study of the filled Julia set of each complex 
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quadratic polynomial reduces to that of an associated polynomial of the type f c . To'thi:, 
end, let p, q, and j denote three complex coefficients, with p ± 0. and let 


/ : C — C, f(u-) = pu' 2 + qw + j. 


(1) Determine complex numbers a, b , and c, in terms of p. q , and j , such thai for 
each u 1 € C 

pu’ 2 + qu.• + j = f(au.’ + b) 2 4- c — f>^j /a. 


(2) Consider the function 


h : C —► C, h(u?) = au' + b. 


Solve h(u') — z for tr to obtain a formula for the inverse function h~ l : C —• C. 
h _1 (z) =? Verify that h(h~ l (z)) = z and h~ l (h(w)) =* w. 

(3) Prove that / = h -1 o f e o h\ thus. f(w) = h~ l (f e (h(w))) for every u’ £ C. 

(4) Deduce that (ir„) = (/ en (u’o)) remains bounded if, but only if, (r n ) = {f° n (h{u-Q })) 
remains bounded. 

(5) Conclude that the filled Julia set of / is the image under h~ l of the filled Julia set 
of f c . Multiplication by a rotates the Julia set by Arg(a) and expands (if jol > 1) 
or contracts (if Jaj < 1) it by a factor of |a|. Adding b translates the Julia set by fc. 

Project 2. 

(1) With c = — 2 the sequence becomes 0, /_ 2 ( 0 ) = 0 2 -2 = -2. 2 (—2) = (—2)*-2 = 

2, /_ 2 ( 2 ) = 2 2 — 2 = 2, and thereafter the sequence stays at 2 forever: in particular, 
it remains bounded, and thus -2 € M. 

(2) Suppose that |cj > 2 and set d = |c| - 2 > 0, so that jc| = 2 + d. Then u'q = 0. 
uq = /c(0) = 0 2 + c = c, u '2 = / c (u'j) = c 2 + c, and 


|u' 2 | = |c 2 +c| = |c(c+1)[ = |c|-|c+l| > |c|-(|c|-l| = |c|-|(2+d)-l| = |c|(l+d) > \c\+2d. 


Suppose that |u> n | > |cj + 2 ■ 4"~ 2 d, which holds for n = 2. Then 


kn + il = |u’ 2 + c| > j|u' n | 2 - |c|! > (|e| + 2 • 4 n-2 d)‘ - \c\ 


= ||c[ 2 + 2|c|2 ■ 4 n ~ 2 d + (2 • 4 n-2 d) 2 - |c|| 

> l c l 2 - |c| + 2 - 2 • 2 • A n ~ 2 d > |cj + 2 ■ 4”" J d. 

Therefore, |ir n | > |c|-f 2-4"~ 2 d, which diverges to infinity. This proves that if |c| > 2 
then c £ M. 

(3) Show that if c £ „Vf, then c € M. 

(4) Since c $ M if jc| > 2, you need only test those complex c for which |cj < 2. 
Second, for a popular algorithm, prove the following proposition. 

If > max{2,(c|} then {u.* n ) diverges to infinity. 
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Proof. Case 1: |c| < 2. If |c| < 2 then max{2, |c|) = 2. Set d — |ir„| - max{2. [c ) = 
|u’ n | - 2 > 0, so that |u.‘„| = 2 + d. Then 

k*+il = |w*S + c\ > ||tr n | 2 - |cl! = (2 + df - lc| = 4 + 2d + d 2 - |c| 

= 2 + 2d + d 2 + (2 - jc|) >2 + 2 d. 

Henceforth, |ir n+ *| diverges to infinity as 2 + 2 k d does. 

Case 2: |c| > 2. If |c| > 2 then max{2, |c|} = |c|. Set d — |tr n | — max{2. c:} = 

|u’n! — Id > 0, so that |u’„l = |c| + d. Then 

k„+iI = k« + d > Ik*! 2 - Id) = (|d + d) 7 - id = Id 2 + 2d|d + d 2 - |cj 

= \c\(\c\ - 1) + 2d|c| + d 2 > \c\ + 4d. 

Thereafter, |u) n +kj diverges to infinity as |c| + 4 k d does. □ 

Third, choose a maximu m number of iterations, N G N. For each complex constant 
c in the closed disc D( 0,2) (or in the square [—2,2] x [—2.2]), compute the first .V 
points of the sequence wo = 0, = 0 2 + c, u>j - c 2 + c, wy = u’ 2 V -i + c • If 

|u\v| > max{2, |c|} then c $ M, by the proposition just proved; color c in white. 

If ju\v| < max{2. jc|} then the results obtained so far do not guarantee that c € M. 
Nevertheless, color c (and c, by symmetry) in black; the accuracy of the resulting picture 
of Mandelbrot’s set depends upon the number of iterations. 

(5) For a better algorithm, prove the preceding proposition with r e r= (1 + ^/l+Tjc i/2 
instead of R c = max{2. |c|). Exhibit S displays working prototypes of programs based 
upon such an algorithm, which produced the pictures in exhibit T. 


Mandelbrot’s set 
on the HP-28C&S 



Tutorial for plotting Mandelbrot’s set 

Mandelbrot’s set on the HP-48SX 

(1) Type, ENTER, and STDre all 
the routines in exhibit 8 

(2) Execute the command Benoit. 

(3) Note that the HP-28 plots the 
upper and lower halves sepa¬ 
rately. 

Exhibit 7. 



For a different program, see William C. Wickes' book [41], pages 337-339. 
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« 1 
WH I l_ 
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REPE 
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IF fi 
THE 

E &fe s 

^ t e 
« CL 
FOR 
— 4" 
SWfl 
FOR 
R 

DU 

-*• 


t e 
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RT 
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LCD 
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Comments on the programs 
for Mandelbrot’s set 

Subroutine tests examines a pixel 
c = (r.y), starting from ro = 
(0,0) and iterating f c ( = ) ~ z 1 + 
c at most twelve times; if ]r t | > 
r c at any iteration, then the sub¬ 
routine leaves c = (r,y) w’hite, 
otherwise it colors c in black. 
ENTER and STOre in 'tests' 

For the HP-48SX only, scene ad* 
justs the screen so that it covers 
the square [ — 2,2] x [—2, 2], which 
contains Mandelbrot's set. 

ENTER and STOre in 'scene' 

For the HP-28Ci:S only, up ad¬ 
justs the screen so that it covers 
the upper half of the square. 

ENTER and STOre in 'up' 

For the HP-28C<k:S only, down ad¬ 
justs the screen so that it covers 
the lower half of the square. 

ENTER and STOre in 'down' 

For the HP-28C<L'S only, sweep 
scans each pixel c = (x,y) in 
D(0,2), which contains Mandel¬ 
brot’s set: sweep computes r c = 
(1+v /rT4i7 j)/2 for each c and 
calls tests. 

ENTER and STOre in 'sweep' 

The main program. Benoit, man¬ 
ages the other subroutines; on the 
HP-48SX, it also performs the tasks 
that sweep does on the HP-28C&:S. 
ENTER and STOre in 'Benoit' 
Tutorial: execute Benoit 
«- on the HP-28C&S, 

on the HP-48SX - 


« 1 * k * STO 

WHILE DUP 
RBS r c x 
k 12 ^ RND 

REPERT 

SQ OVER 
* k 1 1 STO + 

END 

IF RBS rc x 
THEN DUP 

CONJ PIXON 

P I XON 
ELSE DROP 
END » ENTER 

■tests* STO 

« 131 64 x 2 

R+C DUP NEG 

SWRP PDIM # 

131d * 123d 

PD IM .03 125 

* side 1 STO 
»■ ENTER 

* scene ’ STO 

« scene 0 2 

FOR y 

4 y SQ - T 
DUP NEG 

SWRP 

FOR x x y 


R-eC 

DUP 

DUP 

RBS 4 

* 1 

T 1 

-t- 2 

X 

’ r-c 1 

STO 

test 

s side 

STEP 

side 


STEP » ENTER 
‘ Beno 1 1 * STO 


TUTOR IRL 
VRR Benoit 


Exhibit 8. 


















CHAPTER 2 



FRACTALS 

in Real Analysis and Topology 


Summary. The first section presents the construction of von Koch's snowflake and 
shows that it is topologically equivalent to a circle but with infinite length. The second 
section explains the concepts of “Hausdorff dimension" and ‘‘fractal." The present chap¬ 
ter also demonstrates how the description and the explanation of the result of a simple 
construction (by hand or with a computer) require a substantial abstract theory. 

Prerequisites. The construction and the program in the first subsection have no 
formal prerequisites beyond a working knowledge of precalculus. The subsequent explana¬ 
tions of von Koch's snowflake and its fractal nature require the pnor or concurrent study 
of basic topology and real analysis. 
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0. INTRODUCTION 

Benoit B. Mandelbrot coined the word “fractal” in 1975 to name such sets as a circle, 
a segment, all Julia sets (described in the preceding chapter), and von Koch's snowflake 
(pictured in figure 1 and examined in the present chapter). 



Figure 1 . An example of a fractal: von Koch's snowflake. 

Despite its popularity among some mathematicians and scientists, the word "fractal” 
has eluded (and does not need) an accurate definition, perhaps because, though intended 
to designate sets, the word “fractal" refers not to any type of set but to one aspect of all 
sets in the plane and in other spaces: their “Hausdorff dimension" (a concept defined in 
the second section). Thus, studying fractals means studying the Hausdorff dimension of 
any sets. To demonstrate the need for the Hausdorff dimension, the first section constructs 
a particular set, von Koch's snowflake, which has infinite length (a measure of size with 
one dimension) and null area (a measure of size with two dimensions). Consequently, any 
meaningful “measure” of the “size” of such a set must occur with a “fractional dimension.” 
between one and two. 
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1. VON KOCH S SNOWFLAKE 

The present section constructs von Koch's snowflake and investigates its topological 
properties, which reveal that the snowflake is a simple closed curve, topologically equivalent 
to a circle, but with infinite length. 

1.1. Construction and plot of von Koch's snowflake 

The construction of von Koch’s snowflake starts from an oriented segment AB in 
the plane. Informally, an oriented segment AB is a segment of straight line from a point 
.4, designated as the first point, to a point B, designated as the second point. Thus, the 
oriented segment AB. oriented from A to B. differs from the oriented segment BA. 
oriented from B to .4. (Formally, the oriented segment .45 is equivalent to the ordered 
pair (.4.5), which differs from the ordered pair (5,,4) if .4 ^ B.) 


A 


B 


P F 



Figure 2. The basic step of von Koch's construction. 

From the oriented segment AB. the construction of von Koch's snowflake proceeds 
inductively, through a sequence of steps. Basic to the wflole construction, the first step 
transforms an initial oriented segment £ 0 — AB into an oriented polygonal path 
Ei = APQHB. formally equivalent to the sequence of points (A, P, Q. R.B). as shown in 
figure 2, and consisting of four new oriented segments, AP, PQ , QR. and RB. each of 
a length equal to one third that of AB. To specify on which side of the initial segment 
.45 the new vertex Q must lie, adopt the convention that the point Q must be on 
the “right-hand side” of the oriented segment .45; mathematically, this means that the 
ordered basis (AQ,AB) is right-handed, or, equivalently, that the matrix with two rows 
and two columns (AQ, AB) has a positive determinant. Since all four new segments have 
equal lengths, the convention also means that to construct PQ it suffices to rotate .45 
by —7 t/ 3 and then to translate A to P, (Also, applied to the three sides of a triangle to 
produce von Koch’s snowflake, the orientation just described agrees with the orientation 
of the boundary as defined in algebraic topology.) 

After the first step has transformed the initial oriented segment £ 0 = .45 into 
the oriented polygonal path Ei = APQRB , which forms the first stage of von Koch's 
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Exhibit 1 . The first five stages of von Koch's construction, plotted by the HP-2SS 

construction, the second step applies the same transformation to each of the four segments 
AP . PQ< QR< and RB. This produces an oriented polygonal path E? (the second stage 
of the construction) with four times four (sixteen) segments, each with length one third 
that of the preceding ones, in effect one ninth that of AB. and so forth. In general, the k- 
th step produces an oriented polygonal path Ek (the k -th stage) with 4* segments, each 
with length ( i /j) k that of A B, as in exhibit 1. Also, all the 3 x 4 U1 new vertices (new 
corners) produced at the k -th step appear on the same side of the polygonal path as the 
previous ones, thanks to the convention about the orientation adopted in the description 
of the first step. The mathematical “limit’' of von Koch's construction yields one side of 
von Koch's snowflake, denoted by Ex- Exhibit 2 lists a program that continues to execute 
such steps until all the generated vertices exhaust the available memory; exhibit 1 displays 
the consecutive results of such steps. To explain how such a program may work (for any 
type of computer or supercalculator), the following exercises establish some quantitative 
properties of the construction. 


Routine exercises 

Exercise 1. Consider the unit seg¬ 
ment, Eq — [0,1] C R 2 , oriented from 
0 = (0,0) to 1 = (1,0). From that in¬ 
terval, the first step of von Koch's con¬ 
struction produces a polygonal path E\ 

= OPQRl. Calculate the coordinates of 
the vertices P, Q, and R. 

Exercise 2. Instead of the specific segment [0.1], consider any oriented segment 
AB and assume as given the coordinates of its endpoints, A = (x,y) and B = (u.r). 
Calculate the coordinates of the new vertices, P. Q , and H, in terms of j. y. u. 
and u. You may express your results with coordinates, with vectors and matrices, or with 
complex numbers. 
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von Koch’s snowflake 
on the HP-28C&S 

« 'A' STO 'B' STO B A - 
3 / 'C* STO B B C - A C 
turn * - A C ♦ ft I 4 
START line DEPTH ROLLD 
NEXT > 

ENTER 'Koch' STO 


« ClLCD l DEPTH 2 - 
IFERR 

START Koch 
NEXT 

THEN CLEAR 30 SF 
ELSE PRLCD CR CR 
END DEPTH ROLLD 

ENTER ‘pass' STO 


<-13/ INV R*C 1.5 * 
'turn* STO 3 / 2 * DUP 
32 * ‘side' STO INV DUP 
137 * 64 / DUP 0 R+C 
PHAX SNAP R+C NEC PH IK 
(.5.0) (-.5,0) » 

ENTER ' set' STO 


* set CLLCD line 30 CF 
WHILE 30 FC? 

REPEAT pass 
END CLEAR > 

ENTER 'flake' STO 


< DUP2 DUP2 - DUP ABS 
side * CEIL DUP ROT ROT 
/ 'dt' STO 1 SNAP 
START DUP PIXEL dt ♦ 

NEXT DR0P2 » 

ENTER 'line' STO 


TUTORIAL: USER flake 

’ v v v ' 


Comments on the programs 
for von Koch's snowflake 

With points .4 on level 1 and B 
on level 2, this subroutine leaves 
.4, P, Q, R. B on the stack, with 
R = B-(B-.4)/3. 

Q — A — (—1,1/\/3)( B — A)/2. 
P = .4 + (B - A)/3. 

ENTER and STDre in 'Koch' 


With an oriented polygonal path 
(complex numbers) on the stack, 
this subroutine performs one ele¬ 
mentary step of von Koch's con¬ 
struction once for each pair of con¬ 
secutive vertices, thus leaving the 
next entire stage on the stack. 

ENTER and STDre in 'pass' 


This subroutine computes turn 
= (—3, '/3)/2 (thus turn * (.4 - 
B ) = {Q — .4)), the number of pix¬ 
els per unit length side, and the 
dimensions of the screen; it leaves 
the initial stage on the stack. 

ENTER and STDre in 'set' 


This is the main program; it iter¬ 
ates von Koch's construction. 

ENTER and STOre in 'flake' 


Only for the HP-28Ci:S. this sub¬ 
routine draws the segment joining 
two points on the stack. 

ENTER and STOre in 'line' 
Tutorial: execute flake 
in the USER or VAR menu, 


Exhibit 2. 


von Koch’s snowflake 
on the HP-48SX 

< 'A' STO 'B' STO 34* 

3 / 'C STO B 0 C - A C 

turn * - A C + A 1 4 
START 

DUP2 LINE DEPTH ROLLD 
NEXT » 

ENTER ‘Koch' STO 


< ERASE 1 DEPTH 2 - 

IFERR 

START Koch 
NEXT 

THEN CLEAR 30 SF 

ELSE PICT RCL 'snow 1 STO 

END DEPTH ROLLD » 

ENTER ‘ pass' STO 


€-13/ INV R-*C 1.5 * 
'turn' STO 1 3 / INV R*C 
NEC 1 3 / R+C PDIA # 

131d * 131d PDIH (-1,0) 

0 3 / R*C (1.0) (-1.0) > 
ENTER ’set' STO 


< set 30 CF 
HHILE 30 FC? 

REPEAT pass 
END CLEAR » 

ENTER ’flake' STO 


TUTORIAL: VAR flake 
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1.2. Cauchy sequences of continuous functions 

The construction of von Koch's snowflake presented in the preceding subsection pro¬ 
ceeds with a sequence of elementary steps* but the sequence contains infinitely many such 
steps, which means that the final object — von Koch’s snowflake — is a limit. Therefore, 
the analysis of such sets as von Koch's snowflake requires standard facts about limits of 
sequences of functions, which will be reviewed in the present subsection. 

Definition 1. A curve in the plane is a continuous function / : [0. l] —* R : . f — 
/(f) — !’{<)), from the closed unit segment into the plane. (The image /i 0, 1[) C 

R 2 traced by / corresponds to the intuitive concept of curve.) 

EXAMPLE 1. The function s : [0, l] — R 2 with s(f) = (cos(27rf),sin(2::f/) traces 
the unit circle, with center at the origin and radius equal to one. 

EXAMPLE 2, The function /o : [0,1] — R 2 with fo(t) = (t,0) traces the unit 
segment [0.1] x {0} on the first coordinate axis. 


Exercise 3. Write a formula, f\(t) — (A*(f), Y(t )), for the curve OPtJRl corre¬ 
sponding to the first stage of the construction of von Koch’s snow-flake (with four segments). 

Definition 2. For each pair of real numbers a < b the notation C°{[a, fcj. R 2 } 
represents the linear space of all continuous functions from [a, 6] into R 2 . endow-ed with 
the maximum norm, denoted by j| Hoc, and defined for each function / € C 0 ([a.6],R 2 ) 
bv 

11 / 11=0 = max {||/(l)|| : <€[«.,»]}, 

w-ith || || denoting the usual Euclidean norm on R . Thus, ||/||oc represents the greatest 
magnitude that the vector /(f) = (AT(f), y r (f)) reaches as f traces [a, 61. 

REMARK 1. The notation j| ||oc for the maximum norm comes from a standard 

/ t \ 1 / P 

exercise in real analysis, which pertains to the p-norms defined by ||/|j p = ( J a \f\ p ) 

and which shows that ||/||<x> — limp—oc ||/|| p . (See the references by Folland. [22], #7. p 
179, and Rudin, [20], #4, p. 71.) 

Definition 3. A sequence (/„) of functions in C°( [a, 6], R J ) converges to a 
function / € C°([a,6],R 2 ) if, but only if, lim n —« 11/ - /nil* = 0. 

Informally, the definition of “convergence" means that the sequences of vectors (/„(f)) 
= {A" n (/), Vn(f)) converge to /(f) = (X(f),y(f)) at a rate independent of f. 

Definition 4. A Cauchy sequence (/„) of functions in C°([a, 6], R 2 ) is a sequence 
such that for each positive real number e > 0 there exists an integer N t such that if 
n > N t and if m > N t then ||/„ - / m ||oo < f- 

Informally, the functions in a Cauchy sequence get as close as we want to one another, 
but they need not converge to a limit, except in certain particular spaces of functions, for 
instance, in C°([a, 6], R 2 ). 
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Theorem 1. Each Cauchy sequence in C°( [a, 6], R 2 ) converges to a limit. 

Proof. Firstly, observe that for each number t £ a. 6] the sequence of vectors 
(fn(t)) also forms a Cauchy sequence, in R 2 , because ||/„(t) - / m {t)|| < ||/„ - f m \\ < 
if n,m > X e . Consequently, by completeness of the real numbers, each sequence (/„(/!) 
converges to some limit, which we denote by /(*)• That is. f(t) = lim n _^ f n lt). 

Secondly, the limit function / : [a, b] —+ R 2 , /{*) ss limn—^c/nlt) is continuous: if 
r, s £ (a. f>], then the triangle inequality gives 

ll/Cr) - /(J)|| = ||/(r) - Ur) + Ur) - Us) + Us) - !U)\\ 

< l|/(r) - /n(r)|| + ||/„(r) - /„(.)|| + «/.(.) - /(j)||. 

Since (/ n (r)) and (/„($)) converge to /(r) and f(s) respectively, for every t > 0 
there exists an index A' e / 3 such that if n > X t j% then ||/(r) — / n (r)|] < c/3 and 
|[ f{s) - / n (-s)|| < c/3. Moreover, since each /„ is continuous on the closed and bounded, 
hence compact, interval [a. £>], it follows that each /„ is uniformly continuous, which 
means that there exists a positive real number 6 t j 3 such that if r, $ £ fa. 6] and if 
jr - s| < 6 t / 3 then ||/„(r) - )|| < e/3. Therefore, if |r - s\ < 6 r / 3 then 

11/(0 - /(Oil < ||/(r) - /n(r)|| 4 ||/«(0 - / n (j)|| 4 ||/ n (0 - /(s)|| < e/3 4 e/3 4 e/3 = r. 

which means that / is (uniformly) continuous on [a, ft]; thus, / £ C°((a, 6], R 2 ). 

Thirdly, the sequence (/„) converges to its limit / with respect to the maximum 
norm: for each e > 0 and for each t £ [a. 6] there exists an index A T (i */2 such that if 
171 > A'f.r/2 then ||/(f) - /m(OII < c/2, because (/ m (#)) converges to f(t). Also, there 
exists an index X t /2 such that if n, m > A */2 then ||/„ - / m || < c/2, because (/„) is 
a Cauchy sequence. Consequently, if m > X t ,t/2 and if n > A T ( /2 then 

11/(0 -fn(t)\\ < 11/(0 - / m (t)|| 4 ||/ m (t) - /«(0|| < e/2 + c/2 - c. 

Since n does not depend upon t, this inequality means that lim n _oc 11/ - / n |loc = 0. □ 


Exercises 

The following exercises establish the topological properties of von Koch's snowflake. 
To this end. let /„ : [0,1] —► R 1 denote the n-th stage of the construction of von Koch's 
snowflake. For example, /o(0 = (00) represents the initial stage, which traces the unit 
segment, and f\ represents the first stage, with four segments, and so forth. 

Firstly, von Koch's snowflake is a continuous curve. 

Exercise 4. Prove that ||/„+i - /nil* S. l/(3 n 2\/3). Deduce that if 0 < m < n 
then ||fn — fm ||ac < (Vs)" 1 . Conclude that the sequence (/„) converges uniformly to 
a continuous function / : [0,1] —* R 2 , Its image, E 3c = /([0,1]), forms one side of the 
snowflake. 
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Secondly, von Koch's complete snowflake — constructed from each of the three sides 
of an equilateral triangle — is homeomorphic to a circle. 

Exercise 5. By induction, prove that at each stage of the construction the function 
f n is accretive (or expansive), which means that if r.se[a. f>j then ||/„(r) — /„U 1 . > 
|r — sj. Deduce that the limit, /, is also accretive ( ||/(r) —/(-s)|| > |r-s| ), and conclude 
that / is injective: if r ^ s then /(r) ^ f(s). Finally, prove that each side of the 
snowflake is topologically homeomorphic to the unit segment. [0, l], which means that the 
inverse function f~ ] is also continuous on /([0.1]). Thus, the snowflake has topological 
dimension 1. To show that the entire snowflake is homeomorphic to a circle, observe that 
the three sides intersect one another only at their endpoints. 

Thirdly, von Koch's snowflake has infinite length. 

Exercise 6. Prove that the length of the image of /„+j equals 2 * 4 /j that of /„. Then 
find a double sequence of points ((*«,j); = o) * , vou may consider as a sequence 

(fn,;)j = o ea °h stage n € N. such that for each n 

0 — ^n.O ^ ^n.l < ‘ < fn.n —1 ^ ^n,n = 1 

and such that 

lim max {]f nj - • J € {1-,n}} = 0, 

FI — tJC 

n 

lim y:H/U«,i)-/(V,-t)li = =c. 

n — oc ' 

J=1 

This means that for each positive length L > 0 there exists a polygonal path with vertices 
f(t„ j) on the snowflake and with length greater than L. (The curve / — the snowflake 
— is then called non-rectifiable.) 

Lastly, von Koch's snow-flake (which consists only of the “edge") has null area. 

Exercise 7. Prove that the image = /([0, l]) has null area. (Observe that the 
image consists exclusively of the curve /, which traces only the boundary of the flake.) 


2. HAUSDORFF DIMENSION AND FRACTALS 

The existence of sets such as von Koch's snow-flake, which are closed and bounded 
curves covering a null area, like a circle, but with infinite length, unlike a circle, reveals 
the inadequacy of the two standard concepts for measuring the size of planar sets, length 
and area: for von Koch's snow-flake, the one-dimensional measure, length, is infinite, while 
the two-dimensional measure, area, is null. Neither dimension provides a satisfactory 
measure of the size of such a set. This demonstrates the need for a concept of non-integral 
dimensions and for a method of measuring size in such non-integral dimensions, that is. a 
measure between length and area. To develop such a method, the first two subsections, 
below, review the necessary concepts from analysis and topology, and the third subsection 

demonstrates by examples how one might determine the dimension of a set. 
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2.1. The Gamma function and the volume of balls 

Just as any concept of measure requires a unit, the Hau sdorff measure presented here 
relies upon the closed unit balls in R n , denoted by B(0, 1). which consist of all the points 
at Euclidean distance at most one from the origin. The following concepts pertain to the 
n-dimensional volume of the unit ball in R n , which constitutes the unit for measuring the 
size of other sets. 

Definition 5. The gamma function is the function T : ]0,oc[— R defined by 

T(0 = I™ x'-'e^dr 
Jo 


EXAMPLE 3. For example, 

T{1)= f°° x l ~ l e~ z dx = I" t~ x di = —e"X = -0 - (-1) = 1. 

Jo Jo 

From integration by parts follows the recursion formula T(t + 1) = *r(t) for every t > 0: 

T(t + 1) = x t + i ~ 1 c~*dx = x'(—e -r )|{J° — ntx , - 1 (-t-’)dx = 0 + tT(t). 

Jo Jo 

From the recursion formula and T{1) = 1 follow 

r(2) = r(i +1) = ir(i) = i-i = i. 

T(3) = T(2 + 1) = 2T(2) = 2-1 = 2, 

T(4) = T(3 + 1) = 3F(3) = 3-2 = 6. and so forth, 

F(n) = T(n + [n - 1]) = (n - l)r(n - 1) = (n - 1)1 for each positive integer n. 

The calculation of values of T for non-integral arguments involves more substantial 
analysis. For instance, multiples of one-half require the following integral. 

Lemma 1. 

v/ff 
2 


i: 


e"“ du = 


Proof. For details, see Robert Weinstock's note [9], which presents an elementary 
proof along the following line of thought. Consider the function / [0,oc[—» R. w-ith 


/(*) 


- 1 


Notice that /(0) = Arctan(l) = ?r/4, that lim t _oc f(x) = 0, and that 


ru. 


du 
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sU-)= /V-’* 

Jo 

and observe that 

-T = lim {f{x) - /( 0 )> = f f'(x)dx - - f f -=g{-Jx)di 
4 *-* J 0 J 0 y/X 

= -2j~g'{:)g(z)dz = 0-(j~e-'‘ 2 du^ . □ 

REMARK 2. The integral just established is the basis of the “normal'" Gaussian 
distribution in probability and statistics. Another proof proceeds from cartesian to polar 
coordinates in a double integral: 




di 




rdr dB 



—e~* 


9 


ds = 


4 


The alternate proof just sketched is an exercise in many calculus texts, but a rigorous 
proof of the necessary ingredients usually appears only in analysis texts. 

EXAMPLE 4. For t — Vs, setting u = y/x, x = u 2 , and dx = 2 u du gives 


rp/l) = / x l l ^e~ x dx = f u~ 1 e~ u 2udu = 2 f e~ u du = 2~- = %/r. 

Jo Jo Jo 

Hence, the recursion formula gives 

T( Va) = r(l + Vs) = ( l h)T('/ 2 ) = v^r/2, 

T(®/ 2 } - T(*fi + 1 ) = ( 3 /a)r(V 2 ) = 3i/*/4. and so forth. 


Exercise 8 . Calculate T(3.5). 


Definition 6 . For each positive real number d let Qd = r( V*) rf /T(l + d / 2 )- 
EXAMPLE 5. 

o„ = n v 2 )°/r(i + %) = v^°/r(i) = 1/1 = i. 
o, = r(Va)‘/r(i + 'M = -fin-fill) = 2, 
a 2 = r(Vj) 2 /ru + W A = », 

03 = n >/ 2 ) 3 /r(l + W = (v^) 3 /(3v^/4) = 4jt/3. 

Notice that 

Qo = 1 is the number of point(s) in the singleton { 0 }, 

Qi = 2 is the length of [— 1 , 1 ], 

ct 2 — if is the area of the unit disc, and that 

03 = 4tt/ 3 is the volume of the unit ball in R 3 . 





56 


Chapter 2 FRACTALS 


A theorem from analysis also guarantees that aj equals the d-dimensional volume of the 
unit ball in H d for each integral d (see Folland's Real Analysis text. [22], pages 7C-T7). 
For non-integral values of d, the interpretation of qj as the volume of a ball no longer 
holds, but the numbers aj will provide a unit of volume for all dimensions, consistent 
with that for integral dimensions. 

REMARK 3. (SUPERCALCULATORS) Such supercalculators as the HP-2S. the 
HP-48, and their predecessor the HP-15C, have built-in commands (called FACT or !. by 
analogy to factorials of positive integers) to compute the Gamma function, as demonstrated 
in table 1 . 


Table 1 . Computation of the Gamma function with the HP-28(C£:S) and HP-4S 


Keys 

Comments 

Display 

.6 [ red 1 [ *!>£") 1 FACT 1 

(HP-28) Compute TfO.S }. 

.886226925453 

.sIHTHIiprobilu 

(HP-48) Compute r(0.5). 

.886226925453 

1 red 1 n 1 shift] 1 NUM ! | tji \ 2 1 A 1 

Verify with v/ff/ 2 . 

.886226925455 


2.2. Definition of Hausdorff dimension and measure 

The determination of the dimension and size of a set consists essentially of covering 
the set with countably many smaller sets, each with positive diameter at most 6. and 
then adding not necessarily their length, 6, or area, ir( 6 / 2 ) 2 , or volume, *jz^{6j2) 1 . but 
other powers, 6 d . The analysis of the method just described reveals that for each set there 
exists one (exactly one) exponent d > 0 that yields a measure between zero and infinity. 
That exponent will be called the “Hausdorff dimension” of the set. A detailed exposition 
of these concepts requires the definitions in this subsection. 

The following definitions review the foundational concepts related to the construction 
of the real numbers (see Spivalc’s Calculus text, [7], chapter 28). 

Definition 7. (Maximum and minimum) Consider a subset 5 of the real line. 
If there exists a number M € R such that both M € S and s < A/ for every s € 5. 
then M is called the maximum of S, and is denoted by max(S). (If no such number 
A 1 exists, then S has no maximum.) 

Similarly, if there exists a number m € R such that both m £ 5 and m < s for 

every’ $6 5, then m is called the minimum of 5, and is denoted by min(5). (If no 

such number m exists, then S has no minimum.) 

EXAMPLE 6 . The closed unit segment, S = ( 0 , 1 ], has a maximum. A/ = 1 , and a 
minimum, m = 0 . 

Definition 8. (Upper and lower bounds) Let 5 represent a subset of the real 
line. An upper bound for 5 is a number B € R such that s < B for every j £ 5 

Similarly, a lower bound for 5 is a number b 6 R such that i < s for every s € 5. 

(Some sets have no upper or lower bound.) 
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EXAMPLE 7. The open half line S =] — oc.O; has no lower bound but many upper 
bounds, for instance, B — .007. B = 1 . B = \/2. B = e, B — rt, B = 37. anti 
B = 1,000,000 are all upper bounds for ] — oc, 0 [. 

The following theorem from analysis distinguishes the real numbers from the rational 
numbers (for a proof, see Spivak's text, [7], page 556). 

Theorem 2. If a set S C R has an upper bound, then the set V C R of all upper- 
bounds for S has a minimum, min(L’) € U. Similarly, if a set $ C R has a lower 
bound, then the set L C R of all /outer bounds of S has a maximum, maxi L ) € L 

Definition 9, (Supremum and inflmum) If a set 5 C R has an upper bound, 
then the supremum of S is the minimum upper bound for 5, and it is denoted by 
sup(S) or l.u.b.(S) (for /east upper k>und). If not, then we abbreviate the statement 
“5 has no upper bound" by sup(S) = oc. 

Similarly, if 5 has a lower bound, then the infimum of 5 is the maximum lower 
bound for 5. and it is denoted by inf(S) or g.l.b.(S) (for greatest lower iound). The 
notation inf(S) = — oc indicates that 5 has no lower bound. 

The supremum and the infimum differ from the maximum and the minimum by the 
sole additional requirement that the maximum and the minimum of a set must belong to 
that set, whereas the supremum and the infimum need not belong to that set. 

EXAMPLE 8 . The open unit segment, ] 0 ,l[= {i : x € R and 0 < j < 1 ). has a 
supremum, sup’O. 1 [= 1 — min[l, oc[, and an infimum, inf] 0 ,l[= 0 = max] — oc.Q.. Yet 
jO. 1 [ has neither a maximum nor a minimum, because 1 £0,1 [ and 0 £0,1 [. 


Exercise 9. Prove that if S C T C R, then sup(S) < sup(T) and inf(S) > inf(T). 


The following definitions establish the concept of non-integral dimensions. (See also 
Morgan's introductory text [26] and Federer's reference [25].) 

Definition 10. (Diameter) For each non-empty 5 C R n , the diameter of 5 is 
the “extended" real number (finite or infinite, according to the definition 9) denoted by 
diam(S) and defined by 


diam(5) - sup {||x * y || : x, y € 5} , 

By convention, diam(0) = 0 . 

Definition 11. A cover of a subset E C R n is a family C of subsets of R", the 
union of which contains E\ thus, E C U C. A countable cover of E is a cover of £ 
that consists of countably many subsets of R n , of the type C — {S ; : j € N} with each 
Sj C R" and E C U“ 0 5j. (The definition allows some of the sets S } to be empty, so 
that a countable cover may contain either infinitely or finitely many non-empty sets.) 

The smaller the sets 5j in a cover C of a set, the better the cover approximates the 
size of that set. 
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Definition 12. For each £ C R n . for each non-negative rea] number d. and for 
each positive real number 6 , consider all countable covers C = {5j : j € N) of E where 
every S } has diameter at most <5, and define the number 

'Hd.i(E) = inf ja d ^ ^ diam( Sj ) ^ _ £■ q and diam(5 ; ) < <5 for each j € \| . 

which provides a preliminary estimate of the size of the set E in dimension d. 

Observe that if we decrease 6 then we allow fewer covers (only those with smaller di¬ 
ameters), and, consequently, by exercise 9. the infimum increases. Consequently. Tid.fi E) 
tends to a limit or diverges to infinity as 6 decreases. 

Definition 13. (Hausdorff measure) For each set E C R n and for each non¬ 
negative real number d > 0, the d-dimensional Hausdorff measure of E is the 
extended (finite or infinite) real number denoted by H d {E) € R U {oc} and defined by 

H d {E) = lim n d 6 (E). 

f—0 


Exercise 10. From the definition of the Hausdorff measure, prove that if 0 < p < 
q < o c and if Ti p {E) < oo. then 'H,(£‘) = 0. Thus, if the Hausdorff measure of a set is 
a finite number in dimension p, then it equals zero in all higher dimensions q > p. 

Exercise 11. Prove that if 0 < p < g < oc and if H q {E) < oc, then K P {E) = oc. 
Thus, if the Hausdorff measure of a set is a finite number in dimension q . then it is infinite 
tn all lower dimensions, 

The preceding two exercises reveal that if we measure the size of a set in various dimen¬ 
sions, then w>e obtain zero for all higher dimensions and infinity for all lower dimensions, 
and there can exist at most one dimension between both extremes. 

Definition 14. For each set E C R n define the Hausdorff dimension of E as 

the number denoted by h[E) and expressed by 

h(E) — inf {d : d € [0, oo[ and K^(E) < oc}. 


Exercise 12. From the definition of the Hausdorff dimension and from the results of 
exercises 10 and 11, prove the following alternate expressions for the Hausdorff dimension: 

/tf-E) = inf {d : d € [0, oc[ and 'H d (E) =s 0}, 

h(E) = sup {d : d € [0, oc[ and H d [E) > 0}, 
h(E) = sup {d : d € [0,oc[ and Ti d (E) = oc) . 
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Exercise 13. Prove that if 0 < s < oc and if 0 < < x then s = h(E). 

(Such a set E is then called an s-set.) 

The practical significance of the result of the foregoing exercise is that if we determine 
a real number s such that the s-dimensional Hausdorff measure of E is strictly positive 
(but not infinite), then we have also found the Hausdorff dimension s of E. In view of 
exercises 10 and 11 . it also follows that each s-set has a positive finite Hausdorff measure 
only in dimension d = s, because 'HdiE) = oc for d < s and HdiE) = 0 for d > 
Such a result provides a practical method for calculating the Hausdorff dimension of certain 
sets, as demonstrated in the following subsection. 

2.3. Examples of Hausdorff dimensions and measures 


This subsection demonstrates by examples how to calculate the Hausdorff dimension 
and the Hausdorff measure of a few familiar sets in space. 


EXAMPLE 9. (FINITE SET OF POINTS) Let E = consist of k 

distinct points in R n , Further, let r = min {||xi — xj|| : i^j and i,j € {1 . L}}. 

which represents the minimum distance between two distinct points in E\ thus, if Xj ^ 
xj then ||xj — xj|| > r. For each non-negative dimension d and for each positive 

diameter 6 G| 0 . r[, we may cover the set E — {x"\.xT} with k disjoint balls. 

C = {B(xT, 6/2),_B(xT-6/2)}. centered at each point x*j, with radius 6/2. For such a 

cover C, 




/diamB(Xj. 6/2)\ d 

l-2- ) "“‘gU/ 



Consequently. ’Hd.t(E) < 6 d kad/2 d , since "Hd,i{E) represents the infimum of all such 
sums. Therefore. 


0 S = lim Hd.fiE) < (kaj/2 d ) lim 6 d = J W 2< - * if j - “■ 

i-o o 10 if d > 0. 

Hence, by definition of the Hausdorff dimension, E has Hausdorff dimension h(E) = 0 . 

To calculate the 0 -dimensional Hausdorff measure of E , observe that we never need 
more than it sets to cover £, because one set S } for each point xj in E suffices. 
Consequently, since each countable cover C = {5 ; : j € N} D E contains some set S ; 
containing xj, we need only retain k sets Sj,... 5* and we may discard all the others. 
Moreover, 

^diamfSj)^ ^ ^ ^diam(5j) 

jaeQ 

which means that the infimum over such finite subcovers {Si,...,5*} does not exceed 
(and, consequently, by exercise 9, equals) the infimum over all countable covers. However, 
with d = 0 , all such finite covers yield the same measure, because 

V V *- / j — 1 
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Consequently. 

W 0 .s(£) = inf I Q 0 ^ ( diaini Sj) \ . £ c |J 5 jt diam(Sj) < 6 
{ j=i ^ ' j=i 

Finally, 

Ho(£) = lim *{£) — lim k ~ k. 

0 ' d—0 

Therefore, a set with k points has Hausdorff dimension zero and O-dimensional Hausdorff 
measure equal to k. 

Exercise 14. Determine the Hausdorff dimension, d , and the d-dimensional Haus¬ 
dorff measure, of a countable set of distinct points £ = {xj : j € N] C R n . 


| = inf {it} = k. 


REMARK 4. Notice that the Hausdorff dimension of {x"i.x"L} equals 0 regardless 

of the relative positions of the points xT.xT In particular, this means that the 

Hausdorff dimension of a set f C R", like its topological dimension, does not pertain to 
the global position of £ in R” but to an intrinsic and local characteristic of the set 

£ itself. For instance, for each point xj in £ = {xT.xT} the open ball B(xj.r) 

contains only one point from £, namely the point xj, and in that ball the set £ looks 
like one single point, which has dimension 0 . (In contrast, the dimension of the space 
occupied by the set £ relates to the concept of embedding of £ into R n .) 

EXAMPLE 10 . (STRAIGHT LINE SEGMENT) Let £ = [ 0 J] represent a dosed 
and bounded straight line segment of length C in R", from the origin 0 to a point i 
on the first coordinate axis. For each positive diameter 6 > 0 cover the segment ( 0 , C\ 
by rif, = 4 - 1 sets S } , each of diameter at most 6 (recall that the “floor" notation. 

[f/£j, represents the largest integer that does not exceed £/6 ). For instance, cover [ 0 ,f] 
bv |//<fj consecutive intervals (fr + l)<f] of length equal to <5, followed by one interval 
[[f/(fj<5,£] of length at most 6. For such a covering, 

< a,jr < 0i( i(/6i +1) = a „ (p^- 1 +p) 

As the diameter 6 tends to zero, notice that 

(1) W J ([0,fl)=limK<,(([0,<])<o.lim(L^- 1 + = { !f i = ]! 

>-0\- 2 / ( 0 if i>\. 

Consequently, Wd([ 0 ,£]) = 0 for even' d > 1 , and it follows that fi{[0, £]) < 1. 

To prove the converse inequality, consider any cover C = {S ; : j € N} by sets S } 
each of diameter at most 6. Since the intersections I } — S } n [ 0 ,^] C [0,£] lie on the 
interval [0,f], we may replace each set S } by I } and still cover [ 0 ,£]. Further, for each 
positive real c we may replace each set I } by an open inten*al J } l } that 
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contains I } and such that diam(/j) < diam( J } ) < diamf/j) + r/2 J . Moreover, since the 
closed and bounded interval [0, £] is compact in R n , there exists then a finite subcover 

of [ 0 ,£] by finitely many of the open intervals J } , for instance, by {J]. J k } (after a 

rearrangement to avoid double indices). If k = 1 then a single interval J j =]a 1 ,i 1 [ covers 
[ 0 ,£], and i>] — aj > £. Next, proceed by induction: suppose that for some integer m € V 
the inequality — a^) > 6 — 0 = 6 holds for every cover every closed and bounded 

interval [0.6] by k — m open intervals, and consider a cover of [Q.£] by t = m rl 
open intervals ]aj, 6 i [,... 6 m + 1 [ (which need not lie in this order on [0. f]). Also, 

let a, = max{dj : j € (l,...,m + l]} denote the maximum left-hand endpoint among all 
the intervals that contain £, and remove from the cover that interval J, =]a,, 6 ,[. Then 
the remaining m intervals J } , for j ^ i, cover the subinterval [O.a,] C [ 0 .£] because 
[0,Oi]n]a,,f] = 0. Thus, 

rn + 1 m +1 

- Qj) = y ( 6 j - Qj) + ( 6 , - Q.) > (a, - 0 ) + ( 6 , - a,) = 6 , > £. 
j=i /-i 

j*' 

which completes the induction. Consequently, with d = 1 and Qj =2. and for each 
f > 0 , 


( 2 ) 


*.,([ >2±« 


diamC J } ) - s/2 J 


> £ - r. 


Letting t tend to zero gives Hi,$([ 0 ,£]} > £, and then letting S tend to zero yields 
*i([0.£]) > £. With ’Hj([ 0 .£]) < ct ]£/2 — C in ( 1 ), the two inequalities just obtained give 
MM) = 1 , which, by exercise 13, yields the conclusion: 

A closed and bounded segment [0, £] of length £ has Hausdorff dimension d = 
6 ([0 .£]) = 1 and one-dimensional Ha uj dorff measure equal to its length: Kj([ 0 ,£!) = £. 

REMARK 5. More generally than the preceding examples might suggest, the Haus¬ 
dorff measure coincides with the ordinary measure for measurable sets in integral dimen¬ 
sions. Differentiable curves — for example, circles — have Hausdorff dimension 1 , and 
their one-dimensional Hausdorff measure equals their length; similarly, differentiable sur¬ 
faces — for example, spheres — have Hausdorff dimension 2 , and their two-dimensional 
Hausdorff measure equals their area. (See Morgan's text [26], pages 14-17, or Federer's 
[25], Theorem 2.10.35, page 197.) Moreover, as indicated in remark 4, differentiable curves 
have dimension 1 and surfaces have dimension 2 regardless of their position or twisting in 
R n . Thus, the Hausdorff dimension and the Hausdorff measure provide generalizations of 
the intuitive concepts of integral dimensions and of size. 

EXAMPLE 11 . (VON KOCH’S SNOWFLAKE) Cover von Koch's snowflake. E x . 
with the same triangles used to prove that it has null area, in exercise 7. At the k- th 
stage Ek of the construction, 4* triangles S],...S 4 * each of diameter ( 1 /a) t cover 
the snowflake. Consequently, with 6 k = ( 


^<U k 



= a,4‘((>/ J )*/2)‘ = p 
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Observe that 

/ 4 foe if 4/3 d > 1, which means that d< In 4/In 3. 

lim fj = ( 1 if 4/3 d = 1. which means that d *= In 4/In 3, 

^ * [ 0 if 4/3 d < 1, which means that d > In 4/In 3. 

Consequently, Tid(E) = 0 for every d > In4/ln3. from which it follows that h\Ei < 
In4/ln3. Of course, the result just obtained suggests that h(E) = In 4/In 3. but this 
does not yet prove equality, because it merely leads to Hj(B) < oc. not necessarily 
Hi{E) = oc. for every d < In4/In 3. 

The remainder of the proof follows the pattern of the treatment of Cantor's set in 
Falconer’s Fractal Geometry, [34], page 32. To prove that the Hausdorff dimension of 
von Koch's snowflake equals $ = In 4/In 3, it suffices to prove that it has a positive s- 
dimensional Hausdorff measure (see exercise 13). To this end. consider any cover C ~ 
{Sj : j € N} of the snowflake E x . with non-empty sets S } each of diameter 6 } at 
most 6, with 0 < 6 < 1. For each positive real number f. and for each (non-empty) 
set Sj, select any point xj € S } and embed S } into an open ball B } s= B(Xj.rj). 

centered at xj with radius r } = + c. (Such radius r } exceeds the minimum radius. 

6./y/5+e, necessary for the ball B } to enclose the set S } . See the first research project, 
below.) By compactness of the snowflake (the continuous image /([0.1]} of the compact 
unit segment), there exists a cover of E x by finitely many of the balls just constructed. 
Bj,_ B m . For each such ball B r there exists one integer k } such that 

( 1 ) ( l / 3 )‘ J + 1 < diamB ; cp/j)*'. 

Consequently, the ball B } contains at most one vertex of the k } -th stage B* ; . because 
all vertices of that stage lie at least (l/a)* J apart from one another. At each subsequent 
stage Et. with t < k } , the ball B } contains at most 2 x 4 f ~ k > vertices of E(. because 
Bj contains at most two adjacent segments of length ( 1 /a)* J in Et r To relate the 
number of vertices in B } to the diameter of B } . recall that s = In 4 /In 3, so that 3 s = 4: 
consequently, 

( 2 ) 4 '"* J = 4'4 -t ' = 4 , (3 , )" k ' 4 , (3 " < * j + 1)+1 )* = 4 f 3*(3“ (kj +1) )* < 4 , 3’{diamB ; ) # 

by virtue of the first inequality in ( 1 ). Thus, each ball B ; € {Bi. B m } contains at 

most 2 x 4*3*(diamBj )* vertices of Ei. Since the balls cover Boo, a^d since all the 
vertices of the f-th stage Et remain through the subsequent stages and belong to E x . 
it follows that the m balls Bi,...,B m contain the 4* +1 vertices of E(. Therefore, by 
inequality ( 2 ), and with an inequality allowing for duplication (some of the vertices may 
lie in more than one ball), 

m m 

4 * + 1 < ^ number of vertices in B } < 2 x 4 ? 3*(diamBj)V 

Hence, dividing by 2 x 4*3* and rearranging terms gives 

IW ffl tfl 

(3“* + 4“ { ' +1 >)/2 < ^P(diamBj)* = ^(2diam5 ; + 2e) M = V ^(diamS, + f)V 
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Moreover, dividing by 2* twice, letting t tend to zero and f tend to infinity yields 
4~ ( *' t ‘ 1) /2 = 3~»2 ~ j 2~ j /2 < /diamSA < ^ 

j=l ^ 2 ' j = 0 

Therefore, taking the infimum of all such sums for all covers shows that , H J (E 3C } > 
Qj ^-(»+ i >/2 > 0 ; von Koch's snowflake has positive s-dimensional HausdorfF measure 
with s — In4/ln3, which has the necessary consequence that it has Hausdorff dimension 
In 4/In 3 = 1.261859507.... 

Exercise 15. Identify the intersection of von Koch's snowflake and the initial segment 
[0,1] (it coincides with a set famous in analysis and topology). Then calculate its Hausdorff 
dimension. 

2.4. Fractal coastlines versus planimeters 

In suggesting real examples of fractals, most texts liken them to coastlines, implying 
that coastlines may have a non-integral Hausdorff dimension; in particular, this means 
that such coastlines have infinite lengths. 

Exercise 16. Prove that if a set E (which represents a coastline) has a Hausdorff 
dimension h(E) strictly greater than one. then E has infinite length. 

Vet note that a coastline with infinite length would prevent geographers from using 
Amsler's planimeter to measure areas on maps, for the following reason. The planimeter 
applies Green’s Theorem in the plane, which states that for each planar open set fJ C R‘ 
with a simple, closed, and differentiable boundary dQ. and for all differentiable functions 
P.Q . n - R. the following two integrals agree: 

We may think of Green's Theorem as the Fundamental Theorem of Calculus with two 
variables. In particular, for P(x.y) — —y and Q(x,y) = i, Green's Theorem yields the 

Area Formula, 

Area(n > = / J a dx/ ' dv ~\J = i J J a (|f - ^ ) ix A dy 

= ~(£ (-ydx + xdy). 

Based on the Area Formula, Amsler’s planimeter consists of an articulated two-arm instru¬ 
ment, with one end fixed and the other tracing the coastline on a map, and with the pivot 
between both arms recording the line integral in the Area Formula. (For details, consult, 
for example, the reference by Ronald W. Gat ter dam [5].) 

Exercise 17. Suppose that Cl represents the map of a country which is bordered 
by a fractal coastline 3D with Hausdorff dimension strictly between one and two ( 1 < 
ft(3fl) < 2 ), and. by the preceding exercise, infinite length. Explain why a geographer 
may still estimate the area of the country with the planimeter. Would tracing an infinitely 
long coastline not require an infinite amount of time? 
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3. EXPLORATORY TERM PROJECTS 

The present section proposes two projects that may be attempted in parallel with the 
study of this chapter (hence the phrase “term projects”). Since the results involve only 
elementary complex variables and are known to specialists, the projects do not qualify as 
research. However, some of the details are not readily available in many sources, which 
means that you may have to search the literature for some results and derive others yourself 
(hence the adjective “exploratory"), 

3.1. Other fractals of von Koch’s type 

Project 1. Von Koch's construction generalizes 
to other fractals. For instance, exhibit 3 shows how 
to transform an oriented segment AB into an ori¬ 
ented polygonal path APQRSTUVB by inserting 
seven new vertices, which form eight new edges, each 
of length equal to one quarter that of the initial seg¬ 
ment. One of the new' vertices, 5. lies at the mid¬ 
point of AB, with no kink in the path there. (The 
fractal obtained by iteration of the particular construc¬ 
tion shown in exhibit 3 is called a quadric curve of von 
Koch’s type.) 

(1) To understand von Koch's construction better, 
write and test a program to produce consecutive 
stages of von Koch's snowflake on a different type 
of supercalculator or computer. 

(2) Modify the programs for von Koch’s snow r flake to 
experiment with fractals of von Koch's type but 
of your own design. 

(3) Conjecture the value of the Hausdorff dimension 
for the fractals that you create. 

(4) In appropriate cases, generalize the proof of the 
Hausdorff dimension of von Koch's snowdlake (in 
example 11) to your fractals. 

3.2. The smallest ball containing a bounded set 

Project 2. Recall that the calculation of the Hausdorff dimension of von Koch's 
snowflake involved embedding the bounded sets of the cover into balls with radii slightly 
larger than necessary (see example 11, pages 61 - 63). 

For each diameter d > 0 and for each integer n > 2, determine the optimal radius 
R > 0 such that for every set S C R" with diameter d, the set 5 lies in a ball with 
radius at most R. Note that according to the requirements just stated, the radius R 
may depend upon d and n but not upon 5. 






Exhibit 3. 
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4. RESEARCH PROBLEM 


4.1. The Hausdorff dimension of Julia sets 

The research problem suggested here focuses on the Hausdorff dimension and measure 
of quadratic Julia sets, which were presented in the preceding chapter. For the following 
background, see Falconer's Fractal Geometry. [34], pages 20S-218. 

For “small" values of the complex constant c € C, the Julia set J c of the quadratic 
polynomial f c defined by / c (r) = r 2 + c has a Hausdorff dimension approximated by 

<J c = *(r t ) = i + 

Stated more accurately, this means that there exist a positive radius r > 0 and a power 
series of the type a n lc| n , w ith real coefficients a n € R, which converges to the 

Hausdorff dimension, d c = h(J c ) = a "l c i n ' f° r evtfr >' constant c 6 C such that 

\e\ < r. Moreover, ao = 1, Oj = 0, and 02 = 1/(4 In 2). Also, for jc] < r, the d e - 
dimensional Hausdorff measure of the Julia set J c satisfies the inequalities 0 < Hd t (J c ) < 

DC. 

EXAMPLE 12 . If c = 0 then |c| = |0| = 0 < r and d 0 = h(J Q ) = 1 + 0 + 0 + --- = 1. 
Recall that for c = 0 the Julia set Jq of the squaring function is the unit circle, which, 
indeed, has topological dimension and Hausdorff dimension 1; thus, h(J 0 ) = 1 and 
H.Uo) — 2tr (the circumference of the unit circle). 

REMARK 6, If |c[ < 1/4 then the Julia set J c is a simple closed curve (topologically 
homeomorphic to a circle), like von Koch's snowflake. Also as for von Koch's snowflake, if 
0 < jc| < r then 1 < 1 + |c) 3 /(4ln2) + •■• < 2 and the Julia set J c has a non-integral 
Hausdorff dimension between 1 and 2. 

For “large" values of the complex constant c, more accurately, for |c| > (5 + 2\/6)/4. 
the Hausdorff dimension d c of the Julia set J c satisfies the following two inequalities: 

2In2 < < 21n2 

ln(4[|c|+v^R]) " ln(4[| C |-v^R])' 

Problem. Provide a more accurate estimate of the Hausdorff dimension d c of the 
quadratic Julia set J c and of its d c - dimensional Hausdorff measure. 
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5. SOLUTIONS TO ALL THE EXERCISES 


Exercise 1. P = (i/ 3 , 0 ). R - ( 2 /j,0), an d Q = ( 1 / 2 , - V 3 ^3/2). 

Exercise 2. P s= .4+ '/ 3 (B - A), R = A+ 2 / 3 (B~ A), and Q = P + P - A ) = 

P + ( V 2 - -1/3/2)(P - A). 

Exercise 3. 

f (*, 0) if 0 < t < >/ 3 . 

/(<) = { + V»i0) = ((</2)+ '/«.-(< “ 1 / 3 )v / 3/2) if i/ 3 < f < 2 / 3 . 

I (t, 0 ) if 2/3 < f < 1 . 


Exercise 4. At the n-th stage of von Koch's construction, the polygonal path has 
sides of length (Vs)” and differs the most from the next stage at the midpoint of each 
side, by a distance of { 1 /3)" + 1 v/3/2 - l/(3’*2v / 3)- By the formula for the geometric series. 


l!/n ~ /m Hoc S 11 fm + j — /m+j-l Hoc S 

J=1 



m + j-1) 


J U [in/.- 1 1 - ‘ */3 J „ -3 

2^3 h " 2 ^ I" 1 /* " 473 


< 3 -m . 


Consequently, (/„) forms a Cauchy sequence, which converges to a continuous limit, f. 

ExerciseS. Observe that the assertion holds for n = 0, in which case / 0 (f) = (f, 0) 
and ||/o(r) — / 0 {s)|j = ||(r.0) - (s,0)|| = |r — s\. If the assertion holds for some index 
n = k € N, then J|/*+i(r) - /* + ](s)|| > |j/*(r) - /*{s)|| > |r — sj, because the additional 
vertices introduced by fk+i increase distances further. Consequently, by continuity of 
the norm. ||/(r) -/{j)|| = ||Um n _ 0B (/,,(r)-/ n («))|| = lim,,-*, ||/ n (r) - / n (j)|| > |r-«|. 
Thus, if r ^ s then j|/(r) — /(j)|| > |r — s| > 0 and /(r) ^ f(s). Finally, observe 
that / : [0.1] —* /([0,1]) C R 2 is continuous and injective on its compact domain. [0.1]. 
and that its image is a Hausdorff space, which has the necessary logical consequence (a 
standard theorem in topology) that its inverse, / -1 : /{[0, l]) — [0, l], is also continuous. 

Exercise 6. For each stage of von Koch's construction choose for (Vj)j=c a ’l 
n + 1 vertices present at that stage. Since there are then 4" sides each of length { Ya)", 
it follows that |f nj — = f 1 /*)" tends to zero, as required, and that ||/(f n ; )- 

)|| = (*h) n diverges to infinity. 

Exercise 7. Cover the n-th stage of von Koch’s construction, /„, with 4" isoceles 
triangles each of length ( Va)" -1 &Jid height ( l f 3 ) n \fzj2. so that they also cover all the 
subsequent stages, and, consequently, also the limit. The limit — the snowflake itself — 
has am area that does not exceed the total area covered by such triangles, which equals 


4" x 1/2 x (>/3) n-1 x ( l h) n \/S/2 = • • ■ = 3v / 34 n '73 2(n_1) = 3^ Vs)” -1 ■ 
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and which tends to zero. Therefore, the area of the snowflake melts away to zero. 

Exercise 8. I"(3.5) - TfVs} = T ( 5 /2 + 1) = s / 2 r( V 2 ) = V 2 ■ 3^/4 = 15 /ev^* 

Exercise 9. If j£SCl then s € T and. consequently, inf(T) <s< supiT). 
Thus, inf(T) is a lower bound for 5 and sup(T) is an upper bound for 5, and, 
therefore. inf(T) < inf(S) and sup(S) < sup(T). 

Exercise 10. IfO<p<g<oc and ’H P {E) < oc, factor 6 q ~ p out of "Hq.f. which 
gives 


Hqj — inf < 


)~0 x 


/ diamCSj) 


: E C [J and diam(5j) < 6 foreach j > < 

‘ j=o J 


6 9 ~ p inf 


a, —-—— : E C and diam(5j) < 6 foreach j | — t 9 p 'H p j\ ( E i. 


Consequently, as 6 tends to zero. 


K 


,(E) = lim H 
6—0 


fl.* 


(E)< 


lim £ T ' 
0 


p n p _ f = 0 H„(E) = 0 . 


Exercise 11. Proceed as in the previous exercise, but factor 6 P ~ 9 out of ‘H p _t(E) 
and observe that 6 P ~ 9 diverges to infinity as 6 tends to zero, because p — q < 0. 

Exercise 12. Firstly, recall the definition of the Hausdorff dimension: 


h(E) = inffd : d £ [O.ocf and HdiE) < oc}. 

To prove the first alternate expression, consider any d > h(E) and let p denote any 
number such that 0 < p < d < oc, for example, p = ( h(E ) + d) /2. Since h(E) < p the 
definition of h(E) shows that H V (E) < 00 . Then exercise 10, with d — q, shows that 
'Hd(E) ~ 0. which thus holds for even' d > h(E). Yet Hd(E) — x for every d < h(E). 
again by definition of h(E). Therefore, h(E) = inffd : d € [0, oof and Hd(E) = 0}. 

For the second and third expressions, observe again that "Hd{E) = x for every 
d < h(E ) and that Hd(E) = 0 for every d > h(E), by the argument just presented 
for the first expression. Consequently, ME) = supfd : d € [0, cc[ and 'Hd(E) > 0}, and 
fc(E) — sup{d : d £ [ 0 , oof and Hd{E) = 0 }. 

Exercise 13. Suppose that 0 < s < oc and that 0 < 'H t {E) < oc. By exercise 10. it 
follows that HdiE) ~ 0 for every d > s. Also, by exercise 11, it follows that Hd(E) s= oc 
for every d < s. Consequently, s = inf{d : d € [0, ocf and 'Hd(E) = 0} = h(E), by 
definition of the Hausdorff dimension. 


Exercise 14. Proceed as with finitely many points, but cover each point x"j € E 
with a ball S } — B(Xj,6/2 J ) (which may also contain other points jT; € E ). Then 
observe that 



diam{Sj) 

2 





if d > 0, 
if rf = 0. 
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Consequently, h(F) =0 and then 'Hq(E) — ^Jl 0 1 = oo. Therefore, a countable set of 
(distinct) point,! has Eausdorff dimension 0 and infinite O-dimensional Hausdorff measure. 

Exercise 15. From the initial stage of von Koch's construction, Eo “ [0.1], the first 
stage replaces the middle third by the polygonal path PQR. thus leaving the two extreme 
thirds on the unit segment, E\ ~ [0. >/a] U [ 3 /a, 1], Then the second stage in turn replaces 
the middle thirds of both intervals in E\ by polygonal paths, leaving on the unit segment 
the four subintervals Ej = [0, J/s] U [ 2 /s> Vs] CJ [ 2 /a, Vs] U [ B /9-l]i and so forth, w-hich yields 
Cantor’s set. As proved in Falconer's Fractal Geometry, [34J, pages 31-32. Cantor's set 
has Hausdorff dimension d = In 2/ In 3 = 0.630 929 753 57., 

Exercise 16. If E has Hausdorff dimension h(E) > 1 then exercise 12. third 
expression, h(E) — sup{d : d 6 [0,oc[ and Tt^iE) * oc}, shows that if d < h(F) then 
Hd(E) = oc. In particular, since 1 < h{E), it follows that the length of E equals 
7fi(F) = oc. 

Exercise 17. If the coastline dfl is a continuous closed curve, of the type /: [0. l] —* 
R 2 , then a standard theorem of analysis guarantees that for each positive real number r 
there exists a differentiable function / r [0,1] —* R 2 such that / ( (0) = /*(1) (which 
ensures that f e also traces a closed curve) and ||/ - f e \\oc < Therefore, a geographer 
may let the planimeter follow’ the rectified coastline f t and obtain an estimate as close 
as desired to the area of the country. 

5.1. Hints for the term projects 

Project 1. Exhibit 4, on the following page, show’s one way of plotting von Koch's 
snow-flake with Mathematics. 

Returning to Hewlett-Packard supercalculators, modify subroutine Koch so that it 
calculates and inserts the desired points on the stack, for instance B. B — {B - A)/ 4. 
B - {(B — .4)/4}(l,1), B - {{B - A)/4}(2,1), (A + B)/ 2, (.4 + B)j2 + (B - A)(0. l)/4. 
(.4 + B)f 2 + (B - .4)( —1,l)/4, A + (B - A)/4, and .4. Then let Koch iterate the loop 
between START and NEXT eight times instead of four. 

The proof of example 11 generalizes to constructions that do not intersect themselves 
and with an equilateral path at each stage. 

Project 2. In the plane R 2 , if a set 5 has diameter d , then there exists one 
smallest closed disc B containing 5, called the circumscribed disc, and that disc has 
radius r with d/2 < r < dj\/Z\ moreover, there exist planar sets S C R 2 wdth diameter 
d such that r = d/\/3, which means that R — d /\/3 (see Yaglom and Boltyanskn s 
Convex Figures, [27], exercise 3, pages 10 5-106, and exercise 6-1. pages 213-215). 

For R n , prove that d/2 < r < dj ^(n + l)/n, and that there exist set s 5 C R n 
wdth diameter d such that r = d/ y/2 (n + l)/n, which means that R = d/ y/2 (n 4- 1 )/n. 
a result suggested by Kit Hanes. 
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flake:*( (* Mathematics program for von Koch’s snowflake 

u « N[Exp[-Pi*I/3]] ; (* rotates by -1/6 turn *) 

vertices={0,1,Conjugate [u] ,0); (* isoleles triangle *) 

u * 1 + u; {* translates *) 

(* maximum number of iterations ') 
1; (* performs iterations w ) 

(* exhibits the result *) 


stages •= 5; 

Do[ stage, {stages} 
exhibit) 


stage:*( (* "stage" goes 

m • 1; n * 1; (* 

sides * Length[vertices]-1; 
points * 4*sides + X; 
new ■ Table[0, {j,points)]; 

Do[ step, (sides) ]; 
vertices m new) 

step ;* ( (*"step" does 

a * vertices[[n]]; 

b ■ vertices[[n+1]]; 

s - <b-a)/3; 

m * 4*n - 3; 

new ( [m] ] « a; 

new[[a+1]] ■ a + s; 

new[[m+2]] * a + u*s; 

naw[[m+3]] ■ b - s; { 

new([m+4]] ■ b; 

m*m+4;n*n + l) ( 


from one stage to the next stage *) 
next vertex and side to process *) 
(* current number of sides *) 

(* current number of vertices *) 

(* array for new vertices *} 

(* one step for each side *) 

(* retains new vertices *) 

an elementary step on one segment’) 
(* old left-hand endpoint *) 

(* old right-hand endpoint *) 

(* new third of old segment *) 

(* index of next new vertex *) 

(* old first point stays fixed *) 

<* first new vertex at third *) 

{* second new vertex rotated *) 
third new vertex at 2nd third *) 

{* old last point stays fixed *} 
next vertex 4 side to process *) 


exhibit :■( (* "exhibit" displays the latest stage built ’) 

n * Length[vertices]; {* number of vertices *) 

image = {* converts complex numbers to graphics w ) 

Tabla[(Ra[varticas[[i]]],Im[varticas[[i]]]>,(i,n>]; 

ListPlot[image, Axes -> Nona, AspactRatio -> Automatic, 

PlotJoined -> True, PlotStyle -> Thickness[.001] ] ) 

flaka (* execute "flake"; reduce the picture to fit the pace w ) 



Exhibit 4. Program in A/atbemati'ca for plotting von Koch’s snowflake, based upon 
the same algorithm as that for the HP-28 and HP-48 in exhibit 2. 
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accretive, 53 

addition (of complex numbers). 4 
algorithm 
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argument (of a complex number), 3 
attracting (fixed point), 29 
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smallest, 68 
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volume of, 55 

basin of attraction of infinity, 28, 40 
bound, 13, 16, 36, 38-39, 43-44 
bounded (set or complex sequence), 9 
boundary, 23, 63 
Butterfly Effect, 2, 35, 42 
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Cantor's set, 68 
cardiod, 41-42 
Cauchy's Inequality, 25 
Cauchy sequence, 51 
chain rule, 40 
chaos. 2, 35, 42 
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closed 

disc, 22 
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addition. 4 
argument, 3 

arithmetic on the HP-28 and HP-48, 7 

conjugate, 8, 37, 43-44 

derivative, 25 

history, 6 

multiplication, 5 

number, 3 
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plane, 3 

square root, 8 , 38 
continuous complex functions. 21 
convergence 

of Cauchy sequences, 52 
of complex sequences, 10 
of sequences of functions, 51-52 
cover, 57 


curve in the plane. 51 
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dendrite. 2, 34, 42 
dense, 32 

derivative (complex), 25 
diameter, 57 
dimension 

Hausdorff. 58 
topological. 53 
Direct Iteration Method. 14 
disc 

dosed, 22 
open, 22 
unit, 22 

disconnected (J 2 / 5 ), 2, 34, 42 
divergence of complex sequences 9 
Douadys rabbit, 2, 34, 42 
dragon. 2, 34. 42 

E 

expansive, 53 

F 

Feigenbaum {Mitchell J.) ( v (Preface) 
field (of numbers), 6 
filled Julia set, 12 

bound for, 13. 16. 36, 38 
Direct Iteration Method, 14 
interior of, 28 

of any quadratic polynomial 35, 42 
program, 15 
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20 , 34 . 39-40 
fixed point, 16 
attracting, 29 
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parabolic, 29 
repelling, 29 
superattracting. 29 

Fixed Point Inverse Iteration Method, 18 
fractal. 47 

G 

Gamma function, 54 
Gaussian distribution, 55 
Green's Theorem, 63 

H 

Hausdorff 

dimension, 56 

of Julia sets, 64-65 

of von Koch's snowflake. 61-63 
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measure, SB 
I 

imaginary part. 7 
infimum 57 
infinite length. S3 
interior, 28 

invariance of Julia sets. 28, 41 
iteration of complex functions, 11 

J 

Julia set. 26, 35, 42-43 
invariance, 28, 41 
J.* 2, 34 
J 2/ *, 2, 34, 42 
J x/ 4, 19, 29, 33, 35 
Jo, 26, 29, 65 
J-3/4, 34, 40 
J_ 2l 26-27, 34 

K 

von Koch 

construction. 48, 66 
program, 50, 69 
snowflake, 46, 47, 50, 69 

Hausdorff dimension, 61-63 

L 

limit, 21, 52. 58 

of complex sequences, 10 
Lorenz (Edward), 2 
lower bound. 56 

M 

magnitude of a complex number, 3 
Mandelbrot, 47 

set, 35-36, 43^15 
program, 45 

maximum of a set 10, 56 
minimum of a set, 22, 55 
modulus (of a complex number), 3 
multiplication (of complex numbers), 5 

N 

non-attracting (fixed point), 29 
Non-Attracting Fixed-Point Inverse Iteration 
Method, 32 
program, 33 

non-rettifiable (curve), 53 
null area, 53, 66-67 
number field, 6 

O 

open 

disc, 22 


interval, 3 
set, 22, 24-25 

Open Mapping Theorem, 26 

ordered pair, 5 

oriented 

segment, 48 
polygonal path. 48 

P 

parabolic (fixed point), 29* 34 
plammeter, 63 
principal argument, 3 

Q 

quadratic polynomial, 35, 42-43 
quadric curve of von Koch's type, 64, 68 

R 

rabbit (Douady's), 2. 34, 42 

real part, 7 

repelling (fixed point), 29 
Reverse Triangle Inequality, 8* 37 

S 

sequence 

of complex numbers, 9 
of functions. 51 
square root (complex), 8 
f-set, 59 

superattr&cting (fixed point), 29 

supremum, 57 

symmetry 

of filled Julia sets, 16 
of Mandelbrot's set, 35-36, 43-45 

T 

to tend to, 22. 157 
infinity, 9, 26 
topological dimension. 53 
topology, 24 
Triangle Inequality. 4 
Reverse. 8, 37 

U 

unbounded (set or sequence), 9 
unit 

ball. 54 

volume of, 55 
circle. 23 
disc. 22 

upper bound. 56 
V 

volume (of the unit ball), 55 

W 

weather, 2, 42 
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]a. fc[ open interval between o and h , 57 

{a, 6) complex number or pair of real numbers, 3 

{a,&,..>} set consisting of a,6,.,* 

Arg principal argument (of complex numbers), 3 
Qd volume of the d-dimensional unit ball, 55 
C cover (of a set by other sets), 57 
C set of all complex numbers ( C - R 2 ), 3 
c complex number used with / c (r) = z 2 + c t 11 
C°([fl,t], R) space of continuous functions, 51 
/ check mark (after verifications) 

D(r,r) open disc {center z and radius r), 22 

D(z.r) closed disc (center z and radius r), 22 

diam diameter of a set, 57 

6 (Greek “delta*') bound on diameters, 56 

d boundary, 23* 63 

z (Greek “epsilon") tolerance. 51 

/' first derivative of the function /, 25 

f D — C function from D to C, 11, 24-25 

f(z) value of the function / at r, 4 

f ( function with / C (c) - z 2 4- c y 11 

f en n-th iteration of a function /, 11 

"floor" function (greatest smaller integer), 60 
T Gamma function (factorial), 54 
7i d d-dimensional Hausdorff measure, 58 
7i d $ estimate of the sue of a set, 58 
h(E) Hausdorff dimension of a set £\ 58 
i complex number [ i = (0,1), i * * 6 7 = -1 ], 4, 9 
Im imaginary (complex) part of a complex 
number [ Im ( 2 , y) ^ y ], 7 
inf infimum of a set, 57 
oc infinity (shorthand for some limits), 9, 28 
J c Julia set of } ty 26 
h\ filled Julia set of / c , 12 
K* interior of K cy 28 
t length, 60 

lim limit, 9* 10* 21, 25, 51 
M Mandelbrot set, 36, 43-45 
max maximum, 10, 13, 16, 56 
min minimum, 22* 56 
|| || norm, 51 

I] ||ac maximum norm of functions, 51 


N = {0, 1.2,. } set of all natural numbers 9 

Q set of all rational numbers 
r modulus of a complex number. 3: 

(’■ = !=!= \A J + sr ) 

R set of all real numbers 3 
R* set of all non-iero real numbers 
ft+ set of all non-negative real numbers 
R+ set of all positive real numbers 
Re bound ( R c = max{2. |c|) } for Julia sets. 13* 
36 

r € better bound [ r c = (1 4 \/l 4 4|C|)/2 ]. 13, 
36 

Re reaJ part of a complex number [ Re (z . y) ^ z ]* 
7 

sign version of the sign of a real number 8 


sign (y) ~ 



if y > 0. 
if y < 0 


sup supremum of a set, 57 
E (Greek “sigma") summation sign, 53 
0 (‘theta") argument of complex numbers. 3 
\i[ absolute value of the number z, 14 
t *-* /(f) assignment of /(f) to f, 51 
(z t y) ordered pair of real numbers. 3 
{X P(X)} set of all objects X with property 
P[X)> 6, 57 

Z = {... f -2,-l,0 T l'2 _} set of all integers 

7 complex conjugate of 2 , 8* 37 

2 . fixed point, 16 

|jr{ modulus of & complex number, 3 

(z*) sequence of complex numbers, 9 

{zk} singleton (set) consisting of only r*, 31. 55 

Q Greek letter “zeta," 5-6 

[41] reference numbered 41 in the references 

C subset, 9 

U union (of sets), 25, 31, 68 
n intersection (of sets), 25, 32 
\ difference (of sets), 26, 28 
B contains, 33 
€ in, 4 
t not in, 17 

yj complex square root, 8 , 37 
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